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It is well established that many characters in mammals are affected by 
naturally occurring variation in the uterus or, to use a more inclusive term, 
the prenatal maternal environment. Although this factor is regarded as 
environmental as far as the embryo is concerned, it is itself presumably sub- 
ject to both environmental and genetic factors (Russell'). 

The work of Wright,? Wright and Chase,* Green ‘ and others,’ who made 
statistical analyses of variation within inbred strains of mammals, clearly 
shows that the characters which they studied, polydactyly, white spotting, 
number of vertebrae, etc., were affected by environmentally determined 
variation in the prenatal maternal environment. 

The effect of genetically determined variation in the prenatal maternal 
environment on character variation has, however, received little attention, 
perhaps because of a lack of easy methods of evaluating it. Evidence that 
this factor may be important is provided by the increasing number of cases 
in which a difference has been observed between reciprocal F; hybrids ob- 
tained by crossing inbred strains. Russell and Green® have reported one 
such case for number of lumbar vertebrae in mice. Similar data from other 
crosses have been obtained and are in preparation for publication. The 
mere discovery of such a difference does not, however, furnish proof of a 
difference between the maternal environments of the parental strains: 
without further analysis the difference between the hybrids could equally 
well be attributed to the egg cytoplasm or, in the heterogametic sex, to sex 
linked genes. 

Two methods have been developed which can be used to test directly for 
differences between prenatal maternal environments. The first has so far 
been utilized for this purpose on a large scale only by Fekete and Little.’ 
In their work, fertilized ova are removed from the oviduct of one inbred 
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strain of mice and transferred, by means of a fine glass pipette, to the uterus 
of another inbred strain. The resulting offspring thus develop from the 
chromosomes and cytoplasm of the first strain and in the uterine environ- 
ment of the second. 

The second method, developed recently by Robertson,* makes ingenious 
use of transplantation of ovaries within an inbred strain maintained by 
forced heterozygosis. In Robertson’s experiments, ovaries from yellow 
(A*A”) mice were transplanted to agouti (A”A”) mice of the same inbred 
strain. By mating these, and unoperated yellow females, with yellow 
males he was able to compare the development of the lethal homozygous 
A*A” type in the uterus of an agouti with its development in the uterus of a 
yellow. 

The two methods differ in that the transplantation of ovaries brings the 
foster mother’s environment into play at an earlier stage. Whether this is 
important or not is not known, but a combination of the two methods on 
the same material could be used to investigate this point. The methods 
differ also in their practical limitations. The small average number of 
offspring obtained from operations with transferred ova would make the 
collection of the large number of individuals often required for the study 
of quantitative variation somewhat tedious. Ovarian transplantation, 
while more promising in this respect, seemed to be limited to studies of 
single factor genetic differences within inbred strains maintained by forced 
heterozygosis. It was presumably not available for the investigation of 
multiple factor differences between strains because of the failure of tissue 
grafts in foreign hosts. That interstrain transplants of mouse ovaries are 
unsuccessful has been shown by Robertson in a total of 36 operations and 
by Mr. Ralph Kellogg at this laboratory in a total of 23 operations. 

An attempt to find a method by which ovarian transplantation could be 
used to measure the effect of multiple genetic differences in the maternal 
environment led to the experiments described in this paper. It occurred to 
us that if ovaries could be transplanted from an inbred strain to the Fi 
hybrid of that and another strain, a type of transplant that is usually suc- 
cessful, this would provide a method for comparing pure strain and hybrid 
maternal environments. 

For the first experiment, strain dba ovaries were transplanted to dba/ 
C3H F, hybrids which were then mated with dba males. These strains 
were chosen because the dba strain (genotype ddbbaa) differs by three re- 
cessive color factors from the C3H strain (DDBBAA), a fact which, as 
shown in figure 1, would help in distinguishing successful operations from 
those cases in which the hybrid host ovary regenerated. Thus offspring 
from the transplanted ovary would all be dilute brown non-agouti (ddbbaa), 
while only one-eighth of the offspring from a host’s regenerated ovary 
would be of this color. In all operations, except those on females Nos. 162 
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and 163 (see table 1), we followed the essentials of the technique described 
by Robertson.’ The left ovary, capsule and oviduct of each host were re- 
moved. Then, 7 to 15 days later, a second operation was performed in 
which the host’s right ovary was removed from its capsule and the donor 
ovary implanted in its place. In female No. 162 the first operation was 
omitted, one host ovary being left intact, and in No. 163 there was no in- 
terval between the two operations. The latter variation in technique is dis- 
cussed later in this paper. 


DDBBAA x ’ ddbbaa 


C3H STRAIN ‘ | dba STRAIN 


dbo STRAIN & Fi HYBRID @ dbo STRAIN OVA 


| 
1/8 ddbbaa 7/8non-ddbboa ddbbaa 


OFFSPRING FROM REGENERATED OVARY OFFSPRING FROM 
TRANSPLANTED OVARY 


ddbbaa x ES + capbee 
RY 





FIGURE 1 


Of the 22 operated hybrid females, 9 had no offspring and one produced 
a single litter which was eaten shortly after birth. The remaining 12 gave 
the results shown in table 1. Females Nos. 17, 23 and 58 clearly show re- 


TABLE 1 
NUMBER AND TYPE OF OFFSPRING FROM MATINGS OF dba MALES witH dba/C3H. F; 
HYBRID FEMALES INTO WHICH dba OvariES Hap BEEN IMPLANTED 
Serial No. of hybrid host 17 23 25 26 45 55 58 65 162 163 193 197 
ddbbaa offspring wae gage BS Re 2 sr a 1 12 9 8 
Non-ddbtaa offspring Gis Bi re a) RRS RS ae Nee 


generation of the host ovary and failure of the transplanted ovary. Fe- 
males Nos. 25, 26, 55 and 162 also indicate failure of the graft, for the pro- 
portion of ddbbaa offspring produced by each of these females is not signifi- 
cantly different from one-eighth. In females Nos. 45, 163, 193 and 197 the 
transplanted ovary was clearly functional and there was apparently no re- 
generation of the host ovary. The remaining female, No. 65, produced 4 
non-ddbbaa offspring which must have come from a regenerated host ovary, 
but the high proportion of ddbbaa offspring, 18 out of 22, indicates that the 
transplanted ovary was also functioning. 

Robertson® does not mention the possibility of regeneration occurring 
along with a functional graft in his results, but his ratio of 38 yellow to 62 
agouti, where a 1:1 ratio was expected from the transplanted ovary, indi- 
cates that regeneration probably occurred in some of his females and a 
ratio of 38 yellow to 52 agouti from those females which produced at Icast 
one yellow offspring, though not significantly different from a 1:1 ratio, 
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does differ in the direction that favors the possibility of regeneration having 
occurred even in the presence of a functional graft. 

In our work, in any case, subsequent experiments have shown that this 
type of result is not uncommon and apparently cannot be avoided even by 
careful operative technique. This is unfortunate because in the above ex- 
periments, for example, even where the success of the transplant was es- 
tablished by the fact that only ddbbaa animals were born, one cannot be 
certain that any particular ddbbaa offspring was not produced from a re- 
generated host ovary. With hybrids from strains differing by a smaller 
number of marker genes there would be even greater uncertainty, as we 
found in a series of operations made by the junior author and Miss Bar- 
bara Perry.{ In this series C57 black (genotype aa) mice were crossed with 
C3H agouti (AA). C57 black ovaries were transplanted to the F; hybrids 
which were then mated with C57 black males. The expected ratio in the 
offspring from a regnerated hybrid ovary is thus one half agouti (Aa) and 
one half black (aa), while the transplanted ovary would yield only black 
offspring. Of the 23 animals that had young, only one gave reasonably 
clear-cut results. This produced 3 agouti and 18 black offspring, indicating 
a combination of successful transplant and regenerated host ovary. Of the 
remaining 22, 8 gave only black offspring (a total of 26) and 14 gave both 
agouti and black (totaling 56 and 70, respectively). The total of 56 agouti 
to 96 black produced by these 22 females indicates that the ovarian grafts 
were functional in at least some, but the number of offspring was not large 
enough to be sure of this in any particular case. 

It is apparent that the ever-present possibility of regeneration of the host 
ovary occurring along with a successful graft is a serious limitation. The 
difficulty, of course, resides in the fact that the hybrid carries the marker 
genes of both parental strains and will segregate each parental combination 
in at least some of its germ cells. To get around this difficulty a method 
was worked out which makes use of an inbred strain of mice which we have 
called ‘‘stock 129” and which was originally obtained from Prof. L.C. Dunn. 
This strain was inbred with forced heterozygosis to carry both chinch'lla 
(c*) and albino (c*) alleles. It also carries pirlk-eye and light-bellied agouti 
genes, but this can be disregarded in the present discussion. It is now main- 
tained, as shown in figure 2 (a), by brother-sister matings of the heterozy- 
gotes (c“c* X c“c*), thus producing, in addition to the heterozygotes, the 
two homozygous combinations cc“ and c*c*. All three types are pheno- 


typically distinguishable. Two ways in which this strain is being used in 
ovarian transplant experiments are shown in figure 2. In figure 2 (5) the 
F, hybrid is obtained by mating any full colored, CC, strain with the stock 
129 cc type. Ovaries from stock 129 c*c* animals are transplanted to the 
Cc* hybrids which are then mated with stock 129 cc” males. Offspring 


from a regnerated host ovary will thus be either Cc” or cc“, genotypes 
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which are phenotypically distinguishable from the c“c* combination pro- 
duced by a successful ovarian graft. The origin of every individual born 
to the hybrid can, therefore, be determined. Figure 2 (c) shows a similar 
method that can be applied to a cross with any albino strain. The system 
of matings with the strains shown in figure 2 can be worked out in other 


ways, but the two examples given serve to illustrate the principle. 
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FIGURE 2 
(a) System of maintaining stock 129. (b, c) Two ways 
in which this stock can be used in ovarian transplantations 
to hybrid females so that offspring of the grafted ovary can 
be distinguished from offspring of a regenerated host ovary. 


There are other possibilities. Thus cc or c’c’ strains could be used 
just as well as the CC and c°c* strains, and strains carrying other genes 
maintained in forced heterozygosis could be used in place of stock 129. 
As a control for possible effects of the operation itself, transplants can be 
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made both ways between the two homozygous types of the 129 stock or 
whatever stock is used in its place. 

In the above systems, pure strain mice are raised in hybrid hosts by tak- 
ing the transplanted ovary, as well as the male with which its host is 
mated, from the same strain. It is, however, equally possible to use males 
from other strains, thus producing hybrid offspring raised in hybrid uteri. 
One case of this is reported below. Of course the choice of strains is again 
restricted to those carrying suitable marker genes. 

Summar:zing the above, it is clear that, within the limitation of having 
to choose appropriate strains of animals, a method is provided by which ovar- 
ian transplantation can be used for comparing pure strain and hybrid ma- 
ternal environments. The method could also be used in other fields of bi- 
ology wherever, for example, an experiment requires a means of distinguish- 
ing between the functioning of host and graft ovarian tissue. 

In our experiments, ovaries from stock 129 have been transplanted to 
hybrids of this stock and four other stocks: C3H, dba, C57 black and Bagg 
albino. The first three are full colored (CC) strains and were used as shown in 
figure 2 (b), while the Bagg albino was mated as in figure2(c). Functional 
grafts have been obtained in all of these hybrids. In most cases the hy- 
brids were mated to stock 129 males, thus yielding pure stock 129 animals 
raised in four different hybrid maternal environments. Some of the 129/ 
C57 hybrids were mated with Bagg albino males to produce, from the 
transplanted ovary, 129/Bagg hybrids developed in 129/C57 uteri. 

Details of technique, and findings in regard to the effect of various fac- 
tors on the success of the transplantations, will be presented in another 
paper, but it seems desirable to mention here the following departure from 
Robertson’s technique. He removed one host ovary 7 to 10 days before the 
transplantation in order to produce hypertrophy of the remaining ovary and 
of its capsule into which the ovary of the donor was to be placed, but he 
questions whether this is necessary and suggests testing it by removal of 
both host ovaries at the time of transplantation. The success obtained with 
female No. 163 (table 1), which was operated upon in this way, shows that 
hypertrophy of the site of implantation is not an essential factor. In view 
of this finding, most of the experiments with stock 129 have been made by 
removing both ovaries of a host and implanting donor ovaries in their place 
in the same operation. In our work, operations have been more successful 
with this than with Robertson’s method, possibly because two graft ovaries 
per operated female are used instead of one. 

At the present time about 40% of over 100 operations made by the senior 
author on stock 129 hybrids have resulted in successful grafts. Nearly one- 
half of these have shown at least some regeneration of host ovarian tissue. 
The number of offspring obtained from the grafted ovaries has been very 
high in some cases. The highest number recorded for a single female is 80 
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This was a 129/Bagg hybrid containing stock 129 ovaries and mated with a 
stock 129 male. The number of young is higher than has ever been recorded 
at this laboratory for a straight stock 129 mating. (The operation was 
made by Miss Winona Hinkleyt who assisted in the early work with the 129 
stock.) 

A large series of experiments is being conducted to compare the maternal 
environments’ of the 129 stock and some of its hybrids, using number of 
vertebrae as a character. It has already been stated that differences in 
number of lumbar vertebrae obtained in reciprocal hybrids of certain strains- 
indicate that this character may be affected by the maternal environment. 

Summary.—The only way, so far reported, of testing directly for the ef- 
fect of multiple gene, or strain, differences in maternal environments is by 
the somewhat tedious process of transferring fertilized ova from one strain 
to another. The work reported here provides a method of using ovarian 
transplantation, for this purpose, at least as far as a comparison of pure 
strain and hybrid maternal environments is concerned. 

Offspring were obtained from pure strain ovaries transplanted to hybrid 
females. Six different host-donor combinations, involving a total of three 
inbred strains and six hybrids, were tried and all were successful. In many 
cases a successful graft was accompanied by a regenerated hybrid host ov- 
ary. In early experiments it was not possible to distinguish offspring of the 
grafted and regenerated ovaries except on a statistical basis, but mating 
plans were devised using strains with suitable marker genes in such a way 
that it is possible to determine the origin of all the offspring by their pheno- 
type. 

* Aided by a grant from the Rockefeller Foundation and by grants to the Roscoe B. 
Jackson Memorial Laboratory from the National Cancer Institute, the Jane Coffin 
Childs Memorial Fund for Medical Research, the Anna Fuller Fund and the Interna- 
tional Cancer Research Foundation. 
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EXPERIMENTS ON SEXUAL ISOLATION IN DROSOPHILA. V. 

THE EFFECT OF VARYING PROPORTIONS OF DROSOPHILA 

PSEUDOOBSCURA AND DROSOPHILA PERSIMILIS ON THE 
FREQUENCY OF INSEMINATION IN MIXED POPULATIONS 


By HowarpD LEVENE AND TH. DoBzHANSKY* 
CoLUMBIA UNIVERSITY 
Communicated July 19, 1945 


Introduction.—When females of Drosophila pseudoobscura and Drosophila 
persimilis are placed together with males of one of these species, it is usually 
found that a greater proportion of the conspecific than of the alien females 
are inseminated. The nature of the stimuli that lead to this result is still 
obscure, although Mayr and Dobzhansky' have shown that the degree of 
the preference varies depending upon the history of the individual flies 
involved and upon the environment. The experiments to be reported in 
the present article are concerned with the effects of the relative numbers 
of the conspecific and of the alien females available to the males on the 
frequencies of homogamic and heterogamic matings. Perhaps the simplest 
of the many hypotheses that may be constructed about the matings which 
take place in mixed populations is that when a female and a male meet, 
there is a certain definite probability that they will mate, depending on their 
respective species but not on what other males and females may be found 
in the same medium. On the other hand, it is also possible that the flies 
are stimulated by the proximity of conspecific individuals of the opposite 
sex in such a manner that they become more receptive to mating with alien 
as well as with conspecific partners than they would otherwise be. Still 
another possibility is that only homogamic matings take place while 
conspecific individuals of the opposite sex are available, and that hetero- 
gamic matings occur only, or predominantly, when the possibility of homo- 
gamic mating is excluded or remote. The above possibilities are not neces- 
sarily mutually exclusive, and various intermediate situations may be 
encountered. It is also fully conceivable that different pairs of species and 
races might behave quite differently in this respect. 

Methods and Procedure.—The orange strain of D. pseudoobscura and the 
Stony Creek strain of D. persimilis were used in all experiments.! In the 
main body of seven experiments, 10 freshly hatched D. persimilis males 
were placed in each vial together with the following freshly hatched females: 
(1) 10 persimilis, (2) 10 pseudoobscura, (3) 10 persimilis and 10 pseudo- 
obscura, (4) 20 persimilis and 10 pseudoobscura, (5) 40 persimilis and 10 
pseudoobscura, (6) 10 persimilis and 20 pseudoobscura, and (7) 10 persimilis 
and 40 pseudoobscura. In other words, only D. persimilis males were used 
in all vials and their numbers were always ten per vial; the numbers and 
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the species of the females were varied. The vials were kept in an incubator 
at 25'/2°C. for approximately 4 days, whereupon the females were dis- 
sected and the presence or absence of sperm in their seminal receptacles 
was determined by microseopic examination. These experiments were 
done in September—October, 1944, with the assistance of Mr. George 
Streisinger, whose help we wish to acknowledge. The results are sum- 
marized as experiments | to 7 in table 1. 


TABLE 1 


NuMBERS OF INSEMINATED FEMALES OssSERVED WHEN 10 D. persimilis MALES WERE PLACED WITH », D 
persimilis FEMALES AND tq D. pseudoobscura FEMALES 
NO. OF 
CONSPE- NO. OF pe 
EXPT. CIFIC ALIEN ——D. persimilis——~ -———D. pseudoobscura——~ 


NUM- PER PER FER- UNFER- FER- UNPER- ISOLATION ISOLATION 
BER VIAL VIAL TILIZED TILIZED TOTAL TILIZED TILIZED TOTAL INDEX RATIO 
x ez Rar Pez Nez-Fez Nez Faz Naz—Faz Naz bz &z 
eeee * over 
1 10 0 115* 28 143 


80.41* 19.59 


2 0 10 28 125 153 
18.30 81.70 
143.67 34.33 27.34 145.66 
3 10 10 141 7 178 39 «. 184 173. 0.5569 0.2846 
79.21 20.79 22.54 77.46 
149.68 41.33 14.27 78.73 
4 20 10 150 41 191 13 80 93 0.6978 0.1780 
78.53 21.47 13.98 86.02 
204.67 168.33 9.67 80.33 
5 40 10 204 169 37311 79 90 0.6347 0.2235 
54.69 45.31 12.22 87.78 
57.03 44.97 22.99 187.01 
6 10 20 60 , 42 10217 193 210 0.7581 ~—« 0.1376 
58.82 41.18 8.10 91.90 
87.72 41.28 70.17 456.83 
7 10 40 90 39 12964 463 527 0.7035 ~=—«0..1741 
69.77 30.23 12.14. 87.86 
8 0 10 49 80 129 


37.98 62.02 
a a 6 128 67 63 130 
95.31 4.69 51.54 48.46 


* For each experiment, the top line gives the ‘‘expected number’”’ of flies on the basis of H2’ (see text), 
the middle line gives the observed numbers, and the bottom line gives the observed percentages. 


Additional data were secured in April, 1945, with the aid of a slightly 
different method. Males and females of D. persimilis and D. pseudoobscura 
were aged in the absence of individuals of the opposite sex for a week to 
ten days, whereupon 10 persimilis males were placed in vials overnight 
with 10 pseudoobscura females or with 10 pseudoobscura and 10 persimilis 
females. (Experiments 8 and 9 in table 1.) 
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Similar experiments were performed on two geographic strains of D. 
prosalians, but these will be discussed in a separate section. 

When freshly hatched females and males were kept together for 3—4 days 
‘ before dissection, the results obtained in different vials of the same experi- 
ment frequently showed statistically significant heterogeneity, while the 
experiments with the aged flies did not. The heterogeneity is, of course, 
a complicating factor in the analysis of the data. There are several 
probable sources of this heterogeneity. Since dissection and examination 
of the sperm receptacles is a rather laborious operation, females and males 
in some vials were kept together somewhat longer than in others. Some 
flies came from culture bottles in the early and others in late stages of 
hatching. Since the experiments of a given series extended for about two 
months, some variation in the food and other environmental factors may 
have occurred. 

Analysis —Throughout this discussion we will assume that there are 
10 males of a given species or strain present in a standard vial, along with 
n, conspecific? females and n, alien females. Of these f, and f,, respectively, 
will be fertilized. Let NV, equal the sum of the », for a given experiment 
(i.e., NV, is the total number of conspecific females in all the vials of a given 
composition), and let V,, F, and F, be the corresponding sums for mg, f, 
and f,. 

If x, is the probability that any particular conspecific female will be 
inseminated during the course of the experiment, then z, will be some 
function of m,, nm, and of #, the time that the flies are left together. We 
denote this function by 7,(7,, %, t). We will include in ¢ also such factors 
as temperature, age of the flies, food, etc., so that ¢ represents what may be 
called the “physiological time.” Similarly 2,(m,, mq, t) is the probability 
that an alien female will be inseminated. Evidently 7, and x, approach 
zero as the number of females per vial becomes large, since there must be 
some upper limit to the number of females a male can inseminate. Be- 
cause of the variability from vial to vial of the per cent fertilized, as above 
discussed, the data do not give very precise information on the form of 
the functions z, and z,. However, from examination of the results for 
the individual vials it is clear that there is no very great variation in 7, 
over the range of values of m, and m, used, with the possible exception of a 
decrease in w, in the experiment with 40 conspecific and 10 alien females. 
Apparently the decrease in percentage of fertilization does not become 
marked until considerably higher densities of females are reached than were 
used here. Since any analytic function is approximately linear for small 
changes in the variables, we may suppose that 7, depends on some linear 
combination of m, and m, when n, and n, do not vary too widely. We may 
represent this combination as n’ = dn, + en, where d and é¢ are constants. 
If, now, e/d = 0, x, depends only on the number of conspecific females, 
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while if e/d = 1 it depends only on the total number of females. Since 
mw, apparently decreases for n, = 40, mg = 10 and not for m, = 10, m, = 40, 
e/d seems to be a positive fraction less than one. 

Stalker* has proposed an “‘isolation index’ 8 (we — ma)/(we + Oa)- 
An equivalent but somewhat simpler index is y = m,/x,. They may be 
interchanged by the formulas 6 = (1 — y)/(1 + vy) and y = (1 — 8)/ 
(1+ 8). Ifonly homogamic matings occur, 8 = 1 and y = 0; if there is no 
discrimination, 8 = Oand y = 1; and if only heterogamic matings occur, 
68 = —land y = ~. It is somewhat easier to interpret y than £8, as the 
probability of an alien female being inseminated is simply y times the 
probability of a conspecific female being inseminated. We shall call y 
the “isolation ratio.” 

We shall now consider the implications of the hypotheses, discussed 
above, about the manner in which the mating preferences become effective. 
If males do not mate with alien females unless conspecific females are un- 
available, we should expect to find no inseminated aliens until nearly all the 
conspecific females have been inseminated. Since this is not so, we must 
look for some further explanation. It might be that only a small pro- 
portion of the females present are receptive at any one time. This seems 
rather farfetched; it might, however, be tested by seeing what happens 
when only one or two males are present. On the other hand, the hy- 
pothesis might be altered to give a hypothesis which we will call Hl by 
supposing that the probability z, of an alien female being inseminated is 
not zero but has some low value when conspecific females are available 
and has a higher value when only aliens are available. If we also suppose 
that the presence of conspecific females excites males and makes them more 
likely to mate, we will have a situation where 7, has one low value in the 
absence of conspecific females, another low value when conspecific females 
are available, and a higher value when conspecific females have been avail- 
able but are so no longer. The balance between these opposing tendencies 
would then determine whether 7, would be increased or decreased when the 
number of conspecific females m, was increased from zero to, e.g., ten. 
Introducing a second subscript to refer to the experiment number, the 
question is whether m3/mg2 is greater than or less than one. The data 
give’ pqs/Paz = 1.23, but this does not differ significantly from one. This 
comparison was repeated, using aged flies to decrease the heterogeneity as 
explained above. When these vials were tested for heterogeneity they 
gave x” = 19.56 with 28 degrees of freedom, for P = 0.85, so that these 
results are much more reliable. This time pj9/p,3 = 1.36 and this is signifi- 
cant, with x? = 4.81 and P = 0.03. 

Before discussing H1 further we will consider an alternative hypothesis, 
H2. We now suppose that the isolation ratio y = m4(%-, Ma, t)/me(Mey Ma; t) 
is a constant independent of m,, m, and ¢. Evidently this cannot hold for 


II 
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all possible 1,, m, and t; for example, if ¢ is very large and m, and m, small, 
then ,, 7, and hence y will tend to 1; however y may be essentially con- 
stant over a wide range of values. We shall not specify any definite forms 
for the functions 7, and 7, but merely suppose that in the xth experiment 
m, has a fixed value 7,, and 7, has a fixed value 7,,. We temporarily ignore 
the heterogeneity of the data. Hypothesis H2 now reduces to H2’, where 
we assume that 7,/x, has a common value 7, for each vial of experiment x, 
and the hypothesis states that all the y, are equal. We represent the 
common value of the y, by y. 

Evidently y, involves a comparison between two kinds of females and 
hence cannot be computed for experiments 1 and 2, and therefore they 
cannot be used in testing H2’. For the other experiments, under the 
assumption only (i.e., when the y, need not be equal), the best estimates 
for 7,2, etc., are p.;, etc., where p,, is the proportion of conspecific females 
fertilized in all vials of the xth experiment; and the best estimate of 7, 
iS Zr = Par/Per The g, are given in table 1. 

When H2’ is true, the p’s are no longer the best estimates. We may 
now obtain z,, from the relation that under H2’, x,,; = yz, but we still 
require joint estimates of y, m3, ™4, ..., %z- The method of maximum 
likelihood’ leads to the set of six non-linear equations 


7 
3 Eee 
Y 1 — YT, 
— a (1) 
|-a+ = ‘Rae wen Fel] =o (6 = SA. sag ioe 


Ber l =~ Te 1 = hee 





It is not practicable to solve these equations exactly, but a method of suc- 
cessive approximations gives the estimate g = 0.1958 and estimates of the 
Ter Which lead to the ‘‘expected values” that are given in table 1. The 
ordinary x? test is not strictly applicable here because the constraints are 
non-linear. However the likelihood ratio statistic,’ A, is easily calculated 
and, for samples as large as this, —2 log, \ should be distributed approxi- 
mately as x? with 4 degrees of freedom.’ 
To compute —2 log, \ we need the expression 


7 
(—2 log. 10) > (( Fer + Foz) log cz + Fes log Yz + (Nez isthe Fz) log 
x=3 
; (1 ors Tex) + (Naz ria ax) log (1 roe Y2T cx) |. (2) 


(All logarithms are to the base 10 unless otherwise indicated.) Now let 
I, be the value of (2) when we replace y, by y and then replace y, 13, 
Ta) -.-,» Ter by the numerical estimates obtained fom equations (1). Also, 
let L. be the value of (2) when we replace ys, ..., Y7, Tea, .--» ™er, by the 
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sample values given in table 1. Then —2 log, A = L; — L, We obtain 
the result —2 log, \ = 9.0, corresponding to a probability of about 0.061. 
Discussion.—In nature, H2 would mean essentially that there is one 
probability of copulation when a male meets one of his own females and 
another (smaller) probability when he meets an alien female, while under 
H1 we have the added complication that his willingness to mate with the 
alien is reduced by the presence of receptive conspecific females. It should 
be noted that we have spoken throughout as if it were the male who did 
the choosing. This has been done solely for economy of words. Actually 
it is not known whether it is the males or the females or both which exercise 
the discrimination. It may be that the alien females reject the males, and 
that H1 means that the males are more likely to persevere if no conspecific 
females are present. The end result will be the same in any case. The 
important thing is that, under H2, , and w, both decrease at the same 
rate, as m, increases, while, under H1, 7, decreases more rapidly than 7, so 
that y decreases. The data also require the assumption that males be- 
come excited if conspecific females are present, and then are more likely 
to mate with alien females. This additional effect could be superimposed 
on either H1 or H2, but would have to be more extreme to override the effect 
of H1. Since there is no other evidence in favor of H1, it seems that the 
simpler hypothesis: H2 should be accepted. We have seen that, ignoring 
heterogeneity, the probability of the data arising by chance under H2’ is 
about 0.06. The heterogeneity would tend to increase the probability of 
more extreme values of the statistic —2 log, A, and hence to further reduce 
the need of rejecting H2’. On the other hand, the heterogeneity would 
also tend to overshadow any real differences in the y, and hence make the 
test less sensitive. Table 1 shows on its face that the g, do vary over a 
range of 2 to 1, but in an irregular manner rather than in the regular way 
which we would expect under H1. The actual variation in the 7, could be 
even greater than this, but it is not likely that very large systematic changes 
occur. All things considered it seems best to accept the simpler hypothesis 
H2, since it is not contradicted by the data, and indeed seems to fit them 
better than H1. In any event the simple hypothesis that males only mate 
with alien females when their own females are unavailable seems untenable, 
and its modification, H1, is simply a more complicated version of H2. 
Results with Two Strains of D. prosaltans.—Experiments were also carried 
out with the same number of flies as in experiments 1 to 7, but using flies 
of the Chilpancingo (Mexico) strain of D. prosaltans instead of D. persimilis 
and flies of the Bertioga (Brazil) strain of D. prosaltans instead of D. pseudo- 
obscura. As shown by Dobzhansky and Streisinger,* males of both these 
strains exhibit a strong preference for mating with Chilpancingo females. 
Freshly hatched flies were kept together in an incubator at 251/2°C. for 
approximately 3 days, whereupon all the females were dissected and their 
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seminal receptacles examined for sperm. These experiments were per- 
formed in August-September, 1944. The results are summarized as 
experiments 1’ to 7’ in table 2. 


TABLE 2 


D. prosaltans. NUMBERS OF INSEMINATED FEMALES OBSERVED WHEN 10 CHILPANCINGO MALES WERE 
PLACED WITH #, CHILPANCINGO FEMALES AND %q BERTIOGA FEMALES 





NO. OF 
CONSPE- NO. OF 
EXPT. CIFIC ALIEN -——CHILPANCINGO———. ————_- BERTIOGA 
NUM- PER PER FER- UNFER- FER- UNFER- ISOLATION ISOLATION 
BER VIAL VIAL TILIZED TILIZED TOTAL TILIZED TILIZED TOTAL INDEX RATIO 
x Nez Naz Fez Nez-Fez Nez Faz Naz-Faz Naz bz &z 
”. 10 0 57* 33 90 
63.33* 36.67 
2’ 0 10 14 86 100 
14.00 86.00 
3’ 10 10 72 27 99 2 95 97 0.945 0.0284 
72.73 27.27 2.06 97.94 
4’ 20 10 116 58 174 1 86 87 0.966 0.0172 
66.67 33.33 1.15 98.85 
5’ 40 10 156 240 396 1 97 98 0.950 0.0259 
39.39 60.61 1.02 98.98 
6’ 10 20 57 22 79 3 152 155 0.948 0.0268 
72.15 27.85 1.94 98.06 
¥” 10 40 64 32 96 12 381 393 0.912 0.0458 
66 . 67 33.33 3.05 96.95 


* For each experiment, the top line gives the observed number of flies and the bottom line gives the 
observed percentages. 


In these experiments the observed frequencies of heterogamic fertiliza- 
tion are so low that their sampling fluctuations would overshadow any 
reasonable changes in the y,. All that can be said is that the data in 
experiments 3’ through 7’ do not contradict H2’. For such low y values 
large numbers of vials would have to be used to obtain good results. The 
only significant result is the considerable increase of heterogamic matings 
when no conspecific females are present. We have pq3'/Pa2’ = 0.147, with 
x? = 7.43, using Yates’ correction for continuity, and P = 0.006. This 
is probably a real effect, although it might be due to the heterogeneity be- 
tween vials. If real, this effect is more consistent with H1 than with H2. 
It could be at least partly explained under H2 as follows. Suppose z, and 
™_, depend on an “effective population density’”’ with an added female 
affecting’ x, in proportion to her chance of insemination; i.e., if, as sug- 
gested above, we consider z, to depend on a linear combination, n’ = dn, + 
en, = d[n, + (e/d)n,], of n, and m,, then e/d will be equal to y. Now if 
m,(n’, t) decreases rapidly as nm’ increases from 0 to 2 or 4, and thereafter 
decreases less rapidly, then if y is small,.as it is here, 10 y will be small, and 
increasing , from 0 to 10 will reduce a, considerably. There is as yet 
no experimental evidence on the behavior of x, and 2, for small , by which 
this conjecture could be tested. 
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Thus the difference in behavior in the two series of experiments need not 
necessarily imply different causal systems. However, it would not be 
surprising if they were different, since in one case we are dealing with sexual 
isolation between two species and in the other with a one-way sexual 
preference between two strains of a single, different species. 

Summary.—Results obtained by placing D. persimilis males with varying 
proportions of D. persimilis and D. pseudoobscura females admit of the 
hypothesis that the ratio of the probability of a heterogamic mating to the 
probability of a homogamic mating is a fixed constant independent of these 
proportions. However, the possibility of some decrease in this ratio when 
many D. persimilis females are present cannot be rejected. Because of 
small numbers of heterogamic matings, similar experiments with two strains 
of D. prosaltans furnish little evidence on this point. For D. persimilis-D. 
pseudoobscura, heterogamic matings are significantly more frequent in vials 
containing 10 females of each of the two species than in vials containing 
10 alien females only, while for D. prosaltans they are significantly less 
frequent. 


* Experimental data by Th. Dobzhansky, mathematical analysis by H. Levene. 

1 Mayr, E., and Dobzhansky, Th., these PROCEEDINGS, 31, 75-82 (1945). 

2 For simplicity, we shall use the terms ‘‘conspecific’’ and ‘‘alien’’ regardless of whether 
the flies differ as to species (D. persimilis and D. pseudoobscura) or merely as to strain 
(geographic strains of D. prosaltans). 

3 Stalker, H. D., Genetics, 27, 238-257 (1942). 

4 Throughout this paper Greek letters are used for population values (i.e., values 
characteristic of a given experimental set-up) and the corresponding Roman letters are 
used for estimates based on the sample of flies observed. Thus the best estimate of 
we is Poe = Fa/Ne. ‘ 

5 Fisher, R. A., Statistical Methods for Research Workers, London, Sec. 53 ff., 2nd 
(1928) and later editions. 

6 Neyman, J., and Pearson, E. S., Biometrika, 20A, 175-240 and 263-294 (1928). 

7 Wald, A., Trans. Am. Math. Soc., 54, 428-482 (1943). 

8 Dobzhansky, Th., and Streisinger, G., these PROCEEDINGS, 30, 340-345 (1944). 


THE THREE-DIMENSIONAL SHAPES OF BUBBLES IN FOAMS 
By EpwIn B. MATZKE 
DEPARTMENT OF BoTANy, COLUMBIA UNIVERSITY 
Communicated July 23, 1945 
_ The three-dimensional shapes of bubbles in foams have been discussed 
by scientists in diverse fields for nearly three centuries, but apparently 


exact studies have not heretofore been made. 
Even in his Micrographia, published in 1665, Robert Hooke,! the first 
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describer of cells, alluded to “frothy bodies’ and ‘‘congeries of very small 
bubbles,”’ while his contemporary, Nehemiah Grew,? in 1682 drew more 
direct comparisons between plant cells and the ‘‘Froth of Beer or Eggs.” 
Ever since, the similarity between cells in tissues and bubbles in foam has 
been stressed. Thus, Holman, as reported by Child,* made a foam which 
he compared with cells in masses, while Robert‘ published a striking series 
of photographs of groups of bubbles simulating the early cleavage stages 
of the gasteropod Trochus. Berthold’ and Errera® both showed how cells 
and soap films have certain features in common. However, until Lewis’: *: ® 
began his painstaking researches on the three-dimensional shapes of cells 
in tissues, no precise information was available on the shapes of either 
cells or bubbles in foam, even though they had often been compared. 

If biologists have been more interested than other scientists in this topic, 
they are by no means the only ones who have discussed it. Lord Kel- 
vin,’ |! the distinguished mathematician and physicist, who was strongly 
influenced by the scholarly volumes on soap films written by the’ blind 
Belgian physicist Plateau,'? published several essays on ‘‘the division of 
space with minimum partitional area,’’ “homogeneously,” ‘into equal and 
similar parts,” ‘all sameways oriented’’ and under conditions of ‘‘stable 
equilibrium.”” Kelvin solved the problem of division of space into equal 
and similar units with stable angles and each with maximum volume and 
minimum surface by describing a 14-faced figure, the minimal tetrakai- 
decahedron, with 8 undulating hexagonal faces and 6 quadrilateral faces with 
the sides curved. This figure has stable angles. A similar plane-faced 
tetrakaidecahedron, with 8 faces which are regular hexagons and 6 which 
are squares, does not have completely stable angles. Kelvin’s figures have 
been widely accepted, and under ideal conditions with perfect spacing, 
they may well be realized. 

Sir James Dewar'*: |‘: © made extensive studies of soap films and bubbles. 
Separate volumes have been written on the subject by Lawrence’ and by 
Boys,” while Courant and Robbins"™ stress the importance of soap-film 
studies in the solution of certain mathematical problems. 

But the only one who has seriously tried to examine the three-dimen- 
sional shape of bubbles in an actual foam is Desch,! who became inter- 
ested in the problem through its bearing on the form of crystal grains in 
solidifying metals, as in marine propellers and in crucible steel. Desch’s 
data are not extensive, and are somewhat at variance with the results pre- 
sented below. 

This same topic has been cf interest also to students of the alveolar struc- 
ture of protoplasm. Biitschli” concluded that protoplasmic alveoli were 
dodecahedral, while Seifriz?! thinks they may be 12-faceted or 14-faceted. 

It is obvious, therefore, that the three-dimensional shape of bubbles in 
foam has been of interest and of philosophical import to physicists, mathe- 
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maticians, metallurgists, cellular biologists and students of protoplasmic 
structure. It is also evident that no one heretofore has succeeded in look- 
ing into a large mass of soap bubbles and fabulating the number and kinds 
of faces on each. 

In order to accomplish this, at least under one set of conditions, a series 
of studies on the three-dimensional shapes of-bubbles in foams have been 
carried out, using a technique different from that of previous investigators. 
The soap solution found most satisfactory was essentially that of John- 
ston,”? and consisted of the following parts: triethanolamine oleate 7.5 g., 
glycerin 34 g., and distilled water 58.5 g. These substances were shaken 
together thoroughly at intervals for several days, and then allowed to 
stand for at least several weeks. 

Each of the bubbles used was made separately by means of a Yale 
tuberculin syringe graduated to '/,0 cc. If the bubble was to be '/, cc. in 
volume, the plunger was set at '/,19 cc., the tip of the syringe was dipped into 
the soap solution, and the bubble was blown by pressing down the plunger. 
Each bubble was placed separately into a cylindrical dish, having an inside 
diameter of 6 cm. and an inside height of 6.5 cm. The inside of the dish 
was moistened beforehand with the soap solution. A glass plate was used 
as acover. After each bubble was ‘“‘blown,” excess soap solution was re- 
moved by wiping the tip of the syringe on cheesecloth. For most of the 
work the dish was not completely filled, so that the bubbles did not come in 
contact with the cover, even when temperature changes resulted in ex- 
pansion. The dish was rotated from time to time, to keep the solution dis- 
tributed as uniformly as possible and to allow the bubbles to slip and be- 
come readjusted. If the bubbles were to be kept for some time, the dish 
was rotated on its long axis on a clinostat. Although the bubbles would 
last for days, they were always studied as soon as possible after they were 
made, usually on the same day. Even with rotation, some breaking did 
occur, so that the mass of foam cells was not completely static, though a 
state of relatively stable equilibrium was achieved. When bubbles of 
varying sizes were used, the large ones tended to increase in volume, and 
the small ones to decrease, over a period of days. If breaking occurred, 
all the bubbles in the dish were discarded. Although only 1300 bubbles 
were used in the tabulations below, a total of approximately 57,200 bubbles 
were made, each individually, and placed separately into the dish in the 
course of the following experiments. : 

The bubbles, in the dish, were examined under a Spencer binocular dis- 
secting microscope with 1X paired objectives and 6X paired oculars. If the 
films were thin, it was possible to focus through a dish of these bubbles and 
with the dissecting microscope to single out any individual one in a mass 
of approximately 1800, and to record exactly the number and kinds of 
faces. 
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The illustrations shown in figures 1 to 9 were made with a Leitz camera 
lucida attached to the dissecting microscope. The drawings were done 
by plotting the corners of each face and connecting them with straight 
lines. In this way the beautiful, delicate curvatures of the faces and edges, 
by which stability of angles is achieved, are not shown; this method was 
adopted for expediency, since to show these curvatures would be exceed- 
ingly difficult; the numbers and kinds of faces are truly represented. 





FIGURES 1-9 


Camera Lucida drawings of bubbles. Xc 3. The second number under each figure 
represents the total number of faces. Curvatures of faces and edges are not shown in the 
drawings. Quadrilateral faces are designated as Q, pentagonal as P, hexagonal as Hx, 
heptagonal as Hp. The combinations of faces in each are as follows: Fig. 1—11 faces, 
3 Q,6 P, 2 Hx. Fig. 2.—10 faces, 4Q,4 P,2 Hx. Fig. 3.—13 faces, 1 Q, 10 P, 2 Hx. 
Fig. 4.—14 faces, 2 Q, 8 P, 4 Hx. Fig. 5.—12 faces, 12 P. Fig. 6.—10 faces, 3 Q,6P, 
1Hx. Fig. 7.—11 faces, 2Q,8P,1Hx. Fig. 8.—30 faces, 13 P, 16 Hx, 1 Hp. Fig. 
9.—33 faces, 14 P, 17 Hx, 2 Hp. 


In the first set of experiments, bubbles of uniform volume were used. 
The dish was filled, nearly to the top, with bubbles of '/1¢ cc. or 1/2 cc. 
With each filling, bubbles of one volume only were utilized—they were not 
mixed. Since almost identical results were obtained with volumes of 1/10 
cc. and !/2 cc., all the data are combined in table 1. 
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It is immediately evident from table 1 that 400 peripheral bubbles— 
including 200 in contact with the walls of the cylinder and 200 on the 
free upper surface—had an average of 10.99 faces, and that the range was 
from 7 to 14. Bubbles with 11 facets occurred more commonly than any 
other kind. In table 1 it can be seen also that there were more pentagonal 
faces (2338) in these peripheral bubbles than all other kinds combined 
(2057). 


TABLE 1 


NUMBER AND KINDS OF FACES OF PERIPHERAL AND CENTRAL BUBBLES OF UNIFORM 


VOLUME 
NO. OF 7-——NO. OF BUBBLES———. 
FACES PERIPHERAL CENTRAL KINDS OF FACES PERIPHERAL CENTRAL 
7 2 mir Triangular 11 ait 
8 13 ie Quadrilateral 1252 866 
9 20 sie Pentagonal 2338 5503 
10 88 Se. Hexagonal 720 1817 
11 154 2 Heptagonal 73 35° 
12 83 73 Octagonal 1 
13 33 179 Totals 4395 8221 
14 7 218 
15 Uda 106 
16 atts 20 
17 Sis 2 
400 600 


Average 10.99 13.70 


The 400 peripheral bubbles here recorded showed 43 different combina- 
tions of faces. The combination occurring with greatest frequency had 3 
quadrilateral, 6 pentagonal and 2 hexagonal faces (3 Q, 6 P, 2 Hx). It was 
found 67 times among the 400 peripheral bubbles tabulated, and is shown 
in figure 1. Normally the outer face was hexagonal, and there were 6 
lateral and 4 basal contacts. The second most common combination in the 
peripheral bubbles (4 Q, 4 P, 2 Hx) is shown in figure 2; it occurred 47 
times. Although certain of these combinations occurred much more fre- 
quently than others—9 of the 43 were found only once—still no one of these 
may be considered as “the type.” 

In addition to those on the surface, 600 ‘‘central’’ bubbles were examined, 
and the data are recorded in table 1. By “central” are meant those that 
were separated from the top, bottom or sides of the cylinder by at least 
three others. An inspection of the table reveals that the number of faces 
of the central bubbles was higher than that of the peripheral, that those 
with 14 contacts were most numerous, that the range was from 11 to 17, 
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and that the average for 600 was 13.70. It is further obvious that a great 
majority of the faces (5503 of 8221) were pentagonal. 

The 600 central bubbles here tabulated were found in 36 different com- 
binations of faces. The combination found most frequently was a 13- 
faceted foam ceil with 1 quadrilateral, 10 pentagonal and 2 hexagonal con- 
tacts. It occurred 118 times among the 600, and is illustrated in figure 3. 
The second and third most common combinations—1 Q, 10 P, 3 Hx, and 
2 Q, 8 P, 4 Hx (Fig. 4)—were both 14-faceted foam cells, and they were 
found 73 and 64 times, respectively. The fourth most common combina- 
tion of the 36 was the pentagonal dodecahedron (Fig. 5) with 12 faces all 
pentagons; there were 50 of these. But, as in the peripheral bubbles, 
no one of these combinations may be singled out as ‘‘the type,” and al- 
though the average number of contacts is close to 14, Kelvin’s tetrakai- 
decahedra, with 8 hexagonal and 6 quadrilateral faces, are not realized un- 
der the conditions present here. 

In a second set of experiments, bubbles of two different volumes were 
mixed together in known proportions. Small bubbles were made with a 
volume of 0.05 cc., and large bubbles with a volume of 0.4cc. The results 
are summarized in table 2. In experiment 1, table 2, 2 small bubbles were 
made for each large one, so that when the dish was nearly filled and ready 
for examination, it contained twice as many small bubbles as large ones. 
In experiment 2, table 2, 8 small bubbles were put into the dish for each 
large one; in experiment 3, table 2, the ratio was 32 small to 1 large. Fifty 
small bubbles and 50 large ones were then studied in each of these 3 experi- 
ments, and the data recorded. Only ‘central’ small and large bubbles 
were used—that is, none was tabulated unless it was at least the fourth 
one from the top, bottom or sides of the dish. During the course of these 
three experiments, in which 150 small and 150 large bubbles were studied, 
the dish was filled 19 times. 

An analysis of table 2 reveals that when small and large bubbles are as- 
sociated together, the number of contacts of the small ones is less than that 
of bubbles of uniform volume (13.70), while the number of contacts of the 
large ones is increased. Thus in experiment 1, table 2, the average of the 
small bubbles was 9.68, while that of the large ones was 20.42. Further- 
more, increasing the ratio of small to large bubbles brought the number of 
contacts on the small ones closer to 13.70, and as the ratio of small to large 
increased, the number of contacts on the large bubbles became higher, 
within the limits of these experiments. On the basis of the data in table 2 
and of the ratios of small to large bubbles used in the three experiments 
(that is, 2:1, 8:1, and 32:1), the average number of contacts of all the cen- 
tral bubbles (large and small) in the dish at one time can be computed; it 
is found to be 13.26 for experiment 1, 13.62 for experiment 2 and 13.71 for 
experiment 3. Therefore, when bubbles of two different volumes are asso- 
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ciated in varying but known proportions, the number of contacts is altered 
in definite ways. 
TABLE 2 
NUMBER AND KINDS OF FACES OF CENTRAL BUBBLES OF Two DIFFERENT VOLUMES 
ASSOCIATED IN VARYING PROPORTIONS 


NUMBER OF BUBBLES 
NO. OF -———-SMALL BUBBLES- -—————-LARGE BUBBLES——— — 
FACES ExPT. 1 EXPT. 2 ExPT. 3 ExPT. 1 EXPT. 2 ExPT. 3 TOTALS 

















7 1 ‘fs a eS bios du 5% 1 

8 + pier ee a gins iad + 

9 17 1 18 

10 18 3 21 

11 8 14 1 23 

12 2 18 16 36 

13 10 16 26 

14 4 9 13 

15 6 6 

16 1 1 

17 1 2 

18 5 5 

19 6 6 

20 14 14 

21 12 12 

22 10 10 

23 0 se 0 

24 2 4 6 

25 7 7 

26 6 6 

27 10 2 12 

28 8 5 13 

29 5 10 15 

30 9 13 22 

31 1 12 13 

32 7 7 

33 1 1 

50 50 50 50 50 50 300 

Average 9.68 11.90 13.20 20.42 27.34 30.06 
-——SMALL BUBBLES—— -——-LARGE BUBBLES 

KINDS OF FACES ExPT.1 eExPT.2 ExPT. 3 EXPT. 1 ExPT. 2 EXPT. 3 TOTALS 
Quadrilateral 159 118 96 115 57 36 581 
Pentagonal 282 370 412 498 625 624 2811 
Hexagonal 43 ‘101 148 288 557 749 1886 
Heptagonal “tee 6 4 113 118 92 333 
Octagonal ie en iti 6 9 2 17 
Nonagonal 8 ie aay 1 1 at 2 
484 595 660 1021 1367 1503 5630 


From table 2 it can be seen that pentagonal faces occurred more com- 
monly than any other kind, constituting very nearly half the total number. 
Hexagons were also fairly abundant. 
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Varying the volumes and ratios of the bubbles results in an increase 
in the number of combinations of faces. Thus the 150 small bubbles were 
found in 35 different combinations, while the 150 large ones occurred in 105 
combinations. Among the small bubbles the combination found most 
often (16 times) had 10 faces, 3 Q, 6 P, 1 Hx; it is illustrated in figure 6. 
Of the two second most abundant, both of which occurred 12 times, one 
had 11 faces, 2 Q, 8 P, 1 Hx, figure 7, and the other had 13 faces, 1 Q, 10 P, 
2Hx. This latter is also the commonest combination in central bubbles of 
uniform volume (Fig. 3). The large foam cell occurring most frequently 
(7 times) is shown in figure 8; it had 30 faces, 13 P, 16 Hx, 1 Hp. One 
bubble with 33 faces was found—the largest number recorded—14 P, 17 
Hx, 2 Hp (Fig. 9). Here again, no one “‘type”’ is present in either the large 
or the small bubbles. Pentagonal faces occur frequently, hexagonal some- 
what less so, but in no uniform pattern. 

A consideration of the data presented in tables 1 and 2 leads to certain 
definite conclusions: 

1. The peripheral foam cells of a mass of bubbles of uniform volume 
have an average of 11 faces. Lewis** had previously recorded 11-faceted 
epidermal cells. 

2. The central bubbles in a foam, in which all are of approximately 
equal volume, have an average of nearly 14 faces (13.70). But pentagonal 
faces predominate, so that Kelvin’s tetrakaidecahedra with 8 hexagonal and 
6 quadrilateral faces are not realized under these conditions. This does 
not mean that Kelvin was incorrect, for he was interested in the theoretical 
division of space homogeneously into equal and similar units, similarly 
oriented and having stable angles. When a large number of bubbles of 
equal volume are placed into a cylindrical dish, so that they are free to 
glide and adjust themselves, they become arranged in a mass of somewhat 
irregular foam cells, approximating, at least, the equilibrium conditions 
outlined by Plateau, but Kelvin’s ‘‘ideal”’ figure is not achieved. This is at 
variance with the conclusion presented in Sir D’Arcy Thompson’s* schol- 
arly volume. 

3. Comparisons may be made between the data presented above and 
similar data on lead shot of uniform volume (Marvin™) and of varying 
volumes (Matzke*) compressed to eliminate interstices. In such a system, 
starting with spherical shot, the effects of surface energy in subsequent 
changes in shape are obviously of minor importance, in contrast to the con- 
ditions prevailing in soap films. The range in number of contacts, in kinds 
of faces, and especially in number of combinations of faces is more re- 
stricted in the foam cells than in the compressed lead shot. 

4. Detailed comparisons will be published elsewhere between the 
complete data on soap films and cells in plant-and animal tissues, described 
by Lewis, Marvin,” * Higinbotham,” Dodd® and Hulbary.*! It may be 
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noted here, however, that cells in essentially undifferentiated tissues are 
intermediate in many respects between bubbles in foam on the one hand 
and compressed lead shot on the other. And frequently the cells approach 
the bubbles in foam distinctly more closely than they do the compressed 
lead shot. This is interpreted as indicating that surface forces are opera- 
tive when the delicate cell walls are laid down in or near the plant meri- 
stems, that surface forces are of relatively less moment in cellular organiza- 
tion than in soap films, and that they constitute one of the factors, but only 
one, in three-dimensional cell shape determination. 
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GENETIC ANALYSIS OF THE INDUCTION OF TUMORS BY. 

METHYLCHOLANTHRENE. XI. GERMINAL MUTATIONS AND 

OTHER SUDDEN BIOLOGICAL CHANGES FOLLOWING THE SUB- 
CUTANEOUS INJECTION OF METH YLCHOLANTHRENE 


By LEONELL C. STRONG* 
FROM THE DEPARTMENT OF ANATOMY, YALE UNIVERSITY SCHOOL OF MEDICINE 


Communicated July 19, 1945 


The injection into mice, or the application onto the skin surface, of 
methylcholanthrene or any one of a large number of carcinogens is followed 
by a great variety of neoplasms depending upon variations of experimental 
technique. The variation in the types-of cancers obtained with the injec- 
tion of methylcholanthrene has been considerably increased by the use of 
hybrid mice (from the F,-Fjs generations following an outcross) when the 
normally expected local tumors at the site of injection have been partially 
suppressed or completely inhibited through genetic selection toward resis- 
tance to such locally appearing tumors.' 

The actual manner by which a carcinogen initiates normal tissues or 
cells to become cancerous is not clear. It is generally accepted, however, 
that the carcinogen or one of its metabolites acts directly upon cells in 
bringing about the neoplastic change. The observation of Mottram that 
carcinogens cause chromosomes to lag in the equatorial plate during mitosis 
is highly suggestive that nuclear changes may be involved in the origin of 
cancer.?. The recent observation of the author that the untreated descen- 
dants of mice developing cancer of the stomach following the subcutaneous 
injection of methylcholanthrene continue to develop the same type of gas- 
tric cancer spontaneously is evidence that this particular neoplastic condi- 
tion may have been hereditarily established througn at least four genera- 
tions by the effect of methylcholanthrene upon the germ plasm.* Linkage 
tests have demonstrated that susceptibility to gastric carconima is deter- 
mined by the effect of a gene on the “brown tagged’’ chromosome.* 

If it could be demonstrated that germinal changes (mutations) can be in- 
duced by methylcholanthrere, it would be strong presumptive evidence 
that perhaps the action of methylcholanthrene (and the other carcinogens) 
in bringing about cancer is the induction of a somatic mutation. This ex- 
planation is more consistent than the assumption that methylcholanthrene 
has one effect on the germ plasm (induction of mutations) and still another 
totally unrelated effect upon somatic tissue in converting them into cancer. 
Some evidence that mutations affecting susceptibility to tumors induced by 
methylcholanthrene do occur in mice following the injection of the carcino- 
gen has already been published.‘ 

It is the purpose of this paper to present evidence that germinal mutations 
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and other suddenly appearing biological characteristics arise in the de- 
scendants of mice that have received a subcutaneous injection of methy]l- 
cholanthrene and that these sudden changes arise more frequently than 
are to be expected by chance alone. 

Resulis.—One hundred fifty-five F, mice have been produced by a cross 
between C;7 (homozygous for dark eye, black, non-agouti and self) and one 
subline of the NHO descent (proved to be homozygous for dark eye, brown, 
non-agouti and self). The above outcross was made in both directions. 
In each case only the NHO parent was injected with methylcholanthrene. 
Of these, 154 had the expected dark eye, black, non-agouti and self char- 
acters, while one was dark eye, black, non-agouti and spotted (proved to be 
dominant in inheritance). The new variant of spotting involved only the 
ventral surface of the body. Thus a dominant mutation of spotting oc- 
curred once in 155 mice. The author has not had dominant spotting (black 
eye white) in his laboratory for 25 years. It is also possible, although not 
definitely proved, that this new dominant spotting is not the same as the 
old black-eyed white character. 

In all, thirteen color mutations have been obtained in the direct descen- 
dants of mice receiving methylcholanthrene. These are: (1) brown to 
black (dominant change), (2, 3, 4, 5) four pink eye mutants from dark eye 
(recessive inheritance), (6) color to albinism (recessive), (7) spotting in- 
volving only the ventral surface of the mouse discussed above (dominant), 
(8) a second type of spotting involving the nose and ventral surface only 
(dominant), (9) a third type of spotting involving only the belt region of 
the mouse (inheritance uncertain), (10) silvering (precocious graying evi- 
dent upon the appearance of the first hair at 5-7 days of age) (inheritance 
uncertain), (11) permanent wave (inheritance uncertain), (12) piebald to 
self (dominant) and (13) gray belly to non-agouti in a form resembling 
black and tan. The thirteen new variants listed here are all proved to be 
bona fide mutations, based upon the observation of the subsequent appear- 
ance of the new forms in the lineal descendants of the mouse which first 
showed the character. The type of inheritance involved, determined by an 
analysis of the data obtained by an outcross to mice of an unrelated strain, 
has not been made in all cases—hence the use of the expression ‘‘type of in- 
heritance uncertain.” 

Seven of these color mutations are repetitions of characters present in 
the author’s laboratory, while six are new ones never having been present 
previous to the injection of methylcholanthrene. 

The mutation rate is difficult to determine. For example, approximately 
12,000 mice have been injected with 1 mg. of methylcholanthrene. Both 
parents of the NHO descent have been injected with the carcinogen for 
fourteen generations. The restrictions of space and time do not permit the 
keeping of all descendants of the experimental animals. In dealing with 
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hybrid mice as was done in this case, it is necessary to exclude all recessive 
variants that may have been due to segregation. A careful analysis of the 
data indicates that in these experimental animals in one restricted group 
the mutation rate was approximately 1 in 400%. In the controls over a 
period of 27 years, the mutation rate involving coat color changes has been 
approximately 1 in 26,000 mice.* 
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Data on the altered susceptibility to local tumors induced by methylchol- 
anthrene are given in thischart. Data on the original brown non-agouti subline 
are given on the solid line; data for the mutated black non-agouti subline are 
given on the dash line. Time in days is expressed along the base line; percent- 
age incidence of induced local tumors is given on the vertical; that is, the data, 
with advancing time, are cumulative. 


Two sudden changes, other than color changes and alterations affecting - 
cancer susceptibility’ have also arisen in the descendants of methylcholan- 
threne-treated animals. These are. (1) the establishment of precocious 
sexual activity—the production of first litters by 42 days of age in one sub- 
line as compared to 78-85 days, the average in the original ancestral stocks 
which gave rise to them and (2) the establishment of a subline which has 
mice that give rise to extremely large first litters (average 10.2 young) as 
compared to 5.6 young in first litters of mice of the ancestral stocks. 

Two other biological changes have also occurred. Whether they arose 
suddenly is not evident. These are (1) a more rapid growth rate in early 
life and (2) the attainment of adult body weight of 55-60 g. as compared to 
26-28 g. for mice of the ancestral stocks. 
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The derived color mutant sublines either (1) retain the degree of can- 
cer susceptibility (latent period measured from the time of injection to the 
origin of the local tumors) possessed by the ancestral stock mice that gave 
rise to them or (2) suddenly acquire a new susceptibility to induced tumors 
far in excess of the original suscpetibility (Chart 1). The impression is 
gained that many new biological characteristics are being acquired at ap- 
proximately the same time. Thus the conclusion is reached that the effect 
of methylcholanthrene upon the germ plasm at any one time is either (1) 
the induction of a single gene mutation or (2) changes far in excess of what 
are to be expected by a single point mutation. 

Summary.—Thirteen mutations involving coat color have appeared in the 
untreated descendants of mice receiving methylcholanthrene. The first 
mutant obtained occurred in a subline of the NHO strain which had been 
subjected to methylcholanthrene in both parents for twelve generations. 
An increased growth rate, the attainment of a larger body size and weight, 
precocious sexual activity, large first litters and increased susceptibility to 
induced local tumors by methylcholanthrene have also been obtained. These 
suddenly appearing color mutations have occurred either alone or in asso- 
ciation with one or more of the changed biological characteristics enumer- 
ated above. Thus there is evidence that methylcholanthrene has affected 
the germ plasm by bringing about germinal or point mutations and perhaps 
other undetermined effects. It is highly probable, therefore, that methyl- 
cholanthrene may also bring about malignancies in tissues by causing so- 
matic mutations to arise in them. 


* This experiment has been made possible by grants from the Jane Coffin Childs 
Memorial Fund for Medical Research and the Anna Fuller Fund. 

1 Strong, L. C., “‘Genetic Nature of the Constitutional States of Cancer Susceptibility 
and Resistarice in Mice and Men,” Yale Jour. Biol. and Med.,17 289-299 (1944). 
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3 Strong, L. C., ‘‘Genetic Analysis of the Induction of Tumors by Methylcholan- 
threne. IX. Induced and Spontaneous Adenocarcinomas of the Stomach in Mice,”’ 
Jour. Nat. Cancer Inst., 5 339-362 (1945). 

4 Strong, L. C., ““Genetic Analysis of the Induction of Tumors by Methylcholanthrene 
VIII. Two Mutations Arising in Mice Following Injection of Methylcholanthrene,”’ 
Arch. Path., 39 232-236 (1945). 
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DAMPING OF EPIDEMIC WAVES 
By EpwIn B. WILSON AND JANE WORCESTER 
HARVARD SCHOOL OF PuBLIC HEALTH 


Communicated July 12, 1945 


Soper! stated: ‘In the foregoing theory . . . infecting power was sup- 
posed to be instant, at the termination of a certain incubation period, 
and it appears that with such a law in operation, an initial upset of steady 
conditions of prevalence will be followed by epidemic waves that propagate 
themselves indefinitely, without diminution of amplitude.’’ And he went 
on to develop another theory in which infecting power begins at the in- 
stant a person becomes infected and is spread over time but with a dimin- 
ishing incidence; on this assumption he finds that there is damping of the 
epidemic waves. His arguments, be it noted, assume that the distur- 
bance from the steady state may be treated as infinitesimal. 

The first conclusion, namely, that, under the assumption that there is a 
fixed incubation period, there will be no damping, is incorrect, and its in- 
correctness derives from a particular one of the approximations Soper 
makes. He starts out with the equations 


C(t) = (S/m) C(t — 7), C=A — dS/dt (1) 


where C is the case rate, m a constant, S the number of susceptibles and 
A the rate of accession of susceptibles supposed steady, and he remarks 
that since the change in S is usually small in the unit interval (the incuba- 
tion period 7) we may write these equations as 


C(t + 7/2) = (S/m) C(t — 7/2), C=A-— dS/dt. (2) 


This change in effect removes the lag due to the incubation period and 
therefrom is derived the absence of damping in the wave. 

If we use D as a symbol of differentiation and ef as the symbol for 
Taylor’s expansion and if we set S = m + y, with y small so that higher 
powers are to be neglected (1) and (2) become, respectively, 


[-D(1 — e~"?) + Ar/mhy = 0, (1') 
2 sinh = + arhy = 0. (2") 


The second equation has a negative root for D? and hence an undamped 
period in the neighborhood of P = 2r+/mr/A. The first equation does 
not have a pure imaginary solution for D, but a solution with a negative 
real part leading to damped harmonic waves of the approximate form 





y = e 44" cos \/A/mr t 
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with period P = 24+/mz/A (approximately?) and damping per period 
of e—*/2.VAr/m, 

We may give the solution for a more general equation in which the 
generalized form of the law of mass-action is used and not only an incuba- 


tion period 7 but a duration o of infectiveness (both constant) are assumed. 
The equation in S is then* 


p 
4 - (#) = (5) [4 + Sit— 1-2) — St —7)] 


and if S = m + y with y infinitesimal we have 
1 —(r+oe)D 1 


If we let a = 7 + o/2 we have 


D? a? o2 D3 Di— Ap 0. 
E -(4+£)p +(2 +) ee ie], = 
When A p/am is small the solution is approximately‘ 

_ 1pA a . [pA 

= PA (14 Sa) + ma (8) 


and from this the period and damping can be read off. 
For the case of Soper’s distributed infectiousness we should write in our 
notation (with p = 1 in the law of mass-action) 


ee PO ix _ aS\ ¢-n/r 
“ wnat (4 a at 


so that new cases at time ¢ < ¢ are effective at time ¢ only in proportion® 
e*-9/"| The mean value of the duration of infectiveness is then 


So bE + f°. & dt = +. 
Now it so happens that® 





1 
So oO" flOdk = Ty f. 


Hence the equation with the generalized mass law would be 


ds _(s\ D 
4 eh ene 


and for S = m + y, with y infinitesimal, 


D Ap Ap, , AP}, _ 
|>-i3 +4! Ae) y= 0 or [D+ 4? +28] y =0. 
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In this case the value of D may be obtained in finite form as 


p= <2 4 iyi? y! sae 
2m mt 4m 
If compared with (3) taking o = 0 and a = 7+ we see that the waves have 
approximately the same period whether we use a fixed incubation period r 
and an infinitesimal duration of infectiousness or whether we use an in- 
finitesimal incubation period with a period of infectiousness which aver- 
ages to be’r, it being understood that the infectiveness dies out exponen- 
tially. The former case has, however, half the damping of the latter in- 
stead of none at all as Soper states. 

It is of further interest in connection with stepwise calculations to 
note that one may not replace derivatives in differential equations with 
differences without becoming liable to modify the indicated damping. 
For example, in illustrating his theory numerically Soper replaces 


dx _ a as y 
dt 


a : —=a4 
as ar dt as ar ar 


for purposes of calculation by 


xy x 
Ax =a -—-a-—-—, Ay=a-—- ——a-—. 
as ar as ar ar 
If we set x = as(1 + u), y = ar(1 + »v), with u and 2 infinitesimal the 
differential equations, of the infinitesimal epidemic (as found by Soper) 
are 
du dv 


‘<= —q=—-9, TT =a, 


dt dt 


the solution of which leads to the equation for D 








sD.+ 1 Bp ae v 
x ap = 0 = s7D?4+7D+1, 
| 
1 Pee i 1 
D= —- — +4Q— — — % 
2s \~ 4s? 


If, however, we replace D by A = e? — 1 we have 


se? —1) +1 1 
—1 re? — 1) 





= 0 = sre”? — (Qsr — re? +57 —7 +1, 


of which the solution is 
1 1 1 . eT | 
D = 5 tog.(1-2 +4) + tan Y= I 
2 gah gf, 
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If s = 68.2 and r = 2 the values of the pure imaginary parts of the two 
expressions for D will be nearly the same and thus the periods will also be; 
but the real part in the second case will be about —'/,, whereas in the first 
it is —'/.,, i.e., twice as much.’ This explains why Soper found numeri- 
cally about 0.80 instead of about 0.58. 

We have therefore seen that the damping on the “instant hypothesis” 
is determined (approximately) by the logarithmic decrement Ap/4m in- 
stead of zero as stated by Soper and by Ap/2m on the “hypothesis of dis- 
tributed infectiousness.”” We have further seen that the damping is 
modified by replacing derivatives by increments. The rough equations 
given in an earlier paper® for the stepwise calculation of an epidemic with 
recruits, when applied to a variety of hypothetical cases, appear to give no 
damping and, as the theory indicates damping, those equations presum- 
ably have actually eliminated the lag to which damping would be due. 
However, it must be admitted that the phenomenon of recurrent measles 
epidemics gives no clear evidence of any damping. This creates something 
of a difficulty with the theory in respect to the prediction of damping and 
throws some doubt on the reality of periods; it is possible that measles 
simply dies out and then returns and under such a hypothesis there would 
seem to be no reason to expect either definite periods or damping to be ob- 
servable by comparing successive epidemics. 


1 Soper, H. E., ‘The Interpretation of Periodicity in Disease Prevalence,” Jour. Roy. 
Statist. Soc., London, 92, 34-61 (1929). 

2 If a = Ar/m and the root for rD be X + iY we have to solve X — Xe—* cos Y— 
Ye-* sin Yta =0,Y- Ye~* cos Y+Xe~* sin Y=O and the solutions are of the 
fom X = —',a[ltaa+t+aa?+..],¥ =a”? [1+ da+da?+...]. The 
value of a for measles is a small number; if m be 51/2 years of recruits to the susceptible 
population and r be half a month, a = 1/139. 

3 See these PROCEEDINGS, 28, 361-367 (1942) and 31, 24-34, 109-116 (1945). 

4 The value of aD may be written as a power series in pAa/m. If we set 


pe a O4% 7 
m’ 2 


o? 1 o? 1 a? ot 
Ma + Nao toot, 
6 + 24a? 24 + 48a? + 1920a‘ 


AD= in| +(5 M- 27+) | _ pire + (3 N= 5 ra + 1) 


this second approximation will give an estimate of the accuracy of the first approxima- 
tion when definite values of o, r, p, A, mare assumed. Asa = + + o/2 we must have 
o/a varying from 0 to 2 as o/r varies from 0 to ~. Soper really assumed o = 0,7 = a 
fortnights for measles. 

5 The new-case-rate (of becoming infectious) at time #, viz., A — dS/dt, can be re- 
garded as the sum of the products of S by the product 


1 ds * dé ds is ( yf ] 
(4 _~&)\ .-d/s ~~ MmOnitnta. = f-)/t 
m (4 dt ), . site m a] (4 dt ), ' ] dt Jel m ° . 


contributed by the cases which become infectious at the time £ previous to ¢, where the 
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contact rate isr = 1/m. The expression is capable of two interpretations, viz., (a) the 
contact rate remains constant but of the number of cases which become infectious at 
the time £ only the fraction e—"/* remain infectious, which seems to be Soper’s hy- 
pothesis and which accords with the way in which the carrier condition in diphtheria falls 
off with the time, at least after the initial infectiveness (see Stallybrass, C. O., Principles 
of Epidemiology, pp. 302 and 617); (b) the cases which become infectious at time é all 
remain infectious but the contact rate r falls off as e€—/7. In the first interpretation 
7 is the average period of infectiousness, in the second it is the period required for the 
infectiousness to fall off to the fraction 1/e of its initial value; in the first we should 
speak of the diminishing incidence of the infectiousness, in the second of its diminishing 
intensity. 

6 See, for example, Wilson, E. B., Advanced Calculus, p. 451, Exs. 17, 18. 

7 If the root for D is —X + iY, the damping in time is e~ *‘ and if X be divided by 2 
the damping in the same time, being e~ *‘/, is the square root of e~*‘. There is, 
however, a further distinction which must be borne in mind. The differential equa- 
tions are in terms of u = log (C/A). If C/A is near 1 so that the epidemic is infinitesi- 
mal, it is not the case-rate C which is directly subject to damping, but the part of it 
which exceeds A as log(C/A) is practically equalto(C — A)/A. If Cisreally near to A, 
the value of C is practically not damped appreciably. With A = 2200, Cy = 6600, 
Soper’s calculation gave him a second peak of C, = 5168 which is the fraction 0.783 
of the first, and may be taken as a sort of damping factor—very likely the sort an epi- 
demiologist would think most natural. However, to compare with the theory of the 
infinitesimal epidemic one would have to take u = log (6600/2200) = 1.099 and u = 
log (5168/2200) = 0.8544. The ratio of 0.8544 to 1.099 is 0.777 which is very near to 
0.783. The theoretical value for the infinitesimal epidemic is 0.584 but when A=e? —1 
replaces D as an operator the theoretical solution gives, not 0.584, but 0.764 which cer- 
tainly is as near as could be expected to the value 0.777 Soper found. 

Of course, an epidemic with C)/A = 3, uw = 1.099 can hardly be considered infinitesi- 
mal; the period will doubtless have increased somewhat (as Soper actually found in his 
calculation) and the theoretical damping will presumably have been somewhat modified. 
It should be noted that with u = 1.099 we have, if we neglect effects of damping, from 
the equation u, — e“! + e” — uw = 0 between the trough x and crest um, the value u, = 
—1.72; the values +1.10 and —1.72 are so far from symmetrically situated with respect 
to u = 0 as to indicate a major departure from the simple harmonic deviation of the 
infinitesimal oscillation. Indeed it can no longer be expected that the damping as 
estimated by the fraction that one crest is of the preceding one and the fraction that one 
trough is of the preceding one should be the same. A satisfactory approximate solu- 
tion of our equation (3) of a previous paper? 


du 


du 


- = PA pult—1/2)—uit1/2)/9 ( — eultt1/) 
dt\:4-7/2 dt 


t—7r/2 m 








would be required probably to discuss periodicity and damping in the case of finite 
epidemics. 
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TOPOLOGICAL METHODS IN THE THEORY OF FUNCTIONS OF 
A SINGLE COMPLEX VARIABLE. I. DEFORMATION TYPES 
OF LOCALLY. SIMPLE PLANE CURVES* 


By Marston Morse AND MAvRICE HEINS 
INSTITUTE FOR ADVANCED STUDY 
Communicated July 25, 1945 


1. Canonical Curves.—Title I refers to the first of a series of four papers 
by the authors under the general heading of ‘“Topological Methods in the 
Theory of Functions of a Single Complex Variable.” An abstract of Paper 
II will follow this abstract. 

In the study of meromorphic functions one comes most naturally to 
locally simple sensed closed plane curves and deformations of these curves 
through families of sensed curves which are uniformly locally simple. A 
closed curve g is locally simple if for some sufficiently small positive constant 
e each point P of g is an interior point of a simple subarc of g whose end- 
points P, and P; satisfy the condition 


|P,P| >e | P»P| > e. 


A set of closed curves whose points P satisfy the preceding conditions for 
one and the same e will be termed uniformly locally simple. A curve may 
have infinitely many multiple points and be locally simple, and there are 
continuous families of locally simple closed curves which are not uniformly 
locally simple. 

We admit deformations through continuous families of uniformly locally 
simple closed curves. 

A first theorem is that any locally simple sensed closed plane curve g can be 
admissibly deformed into a curve C” obtained by tracing a sensed circle 

x = cos 0 y = sin 0 (1.1) 

n times,n = +1, +2, ... or intoa figure eight, C°. Moreover g can be de- 
formed into but one «f these canonical curves. 

In particular the figure eight, C°, traced m times can be admissibly de- 
formed into C°. 

If a locally simple curve is given as the locally 1-1 continuous image of 
the circle (1.1) in the form 


w = f(@) (w = u + iv) 
then as @ increases from 0 to 2x any continuous branch of the function 
are [f(@ + e) — f(@)] (0 <eS 4) 


will change by an integral multiple 2px of 27, independent of e provided e; 
is a sufficiently small positive constant. The integer ? is called the angular 
order of g. : ; 
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It is apparent that the angular order of an admissible curve g is invariant 
under any admissible deformation of g. It is shown that an admissible curve 
of angular order p can be admissibly deformed into the above canonical curve 
C?, including the figure eight C° when p = 0. 

A polygon z is locally simple if no edge has more than one end-point in 
common with the succeeding edge of the polygon. A continuous variation 
of a circular sequence of vertices representing a polygon will define an 
admissible deformation of 7 if successive vertices are bounded from each 
other and if successive edges have at most an end-point in common. We 
admit the addition or removal of a finite number of vertices on edges of a 
polygon at a finite number of moments ¢ of a deformation. 

The first approach to the preceding theorems is by way of lemmas by 
virtue of which any locally simple sensed curve can be admissibly deformed 
into an admissible polygon 7. It is then shown that 7 can be admissibly 
polygonally deformed either into a triangle traced m times, n = +1, +2, 
..., or into two triangles with a vertex in common forming a figure eight. 
The proof is by induction with respect to the number of multiple points. 

2. Products of Deformation Types.—The class of admissible curves 
admissibly deformable into an admissible curve g will be denoted by (g) 
and termed a deformation type 7. Each deformation type includes a 
sensed regular curve. With the aid of such regular curves the product 


T,T2 = (g1) (ge) 


will be defined. Let i, and hz be regular sensed curves in (g;) and (ge), 
respectively, which are positively tangent at a point P. Let Mhz be a 
closed curve obtained by cutting /,; and hf, at P, to form arcs h’ and h’, 
respectively, tracing h’, then h”, and joining the final end-point of h” to 
the initial end-point of h’. We define 7,T2 as the deformation type (/he) 
and show that this type depends only on 7; and 72, and not on P or the 
choice of /, and hz in the types 7; and 7>. 

It is shown that the product T;T> 1s commutative and associative. 

For n positive or negative let g” denote g traced m times, For n = 0 
let g” denote a curve of the type of (g)(g~"). With these definitions it is 
shown that 


ie”) = o) 


where r and m are arbitrary integers. It must be clearly understood that 
the operation ( )” isa priori utterly different from the operation of tracing 
gm times. “tos 

It follows that the deformation types form an Abelian group with respect 
to the product operation, with the type of the figure eight the idem-factor. This 
group 1s tsomorphic with the additive group of integers. 
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3. The Order q.—The order q(w) of a sensed closed curve g with respect 
to point w not on g is well defined. 

If g is locally simple, 1t is shown that the order q(w) possesses an absolute 
bound independent of points w not on g. 

An admissible deformation none of whose curves intersects a fixed point 
O will be called an O-deformation. It will be convenient to take O as the 
origin in the w-plane. 

We need canonical curves with respect to O-deformations, given p and 
q. To that end let C be the sensed unit circle previously defined. Let C, 
be a similarly sensed circle of radius 1/4, tangent to C at (u, v) = (1, 0). 
Let 


CC (q #0,r # 0) (1.2) 


denote the curve obtained by first cutting C and C, at (1, 0), starting with 
the point (1, 0), then tracing C q times in the sense indicated by the sign 
of q, next tracing C; similarly r times, closing the curve at (1, 0). The 
curve C, shall be taken interiorly or exteriorly tangent to C at (1, 0) accord- 
ing as q and 7 have the same or different signs. The resulting curve will 
be locally simple and regular. It is observed that the curve (1.2) has the 
order g with respect to w = 0 and the angular order p. 

There exists an admissible curve g with prescribed order g and angular 
order p. A principal theorem follows: 

Let g be an admissible curve of order q relative to the origin and of angular 
order p. Let the following cases be distinguished 


CaseI g #0 pPp-q#OD 
Case II g #0 P-q=0 
Case III g = 0. 


Then g admits an O-deformation into a curve C*C,?—* in Case I, a curve of the 
form C? in Case II, and a curve of the form C,? in Case III. No one of these 
canonical curves with the orders (p, q) admits an O-deformation into a canonical 
curve with a different pair of orders. 

The relation of this deformation theory to the theory of meromorphic 
functions will become clear in Paper II. 


* Morse, M., and Heins, M., ‘‘Topological Methods in the Theory of Functions of a 
Single Complex Variable: II. Boundary Values and Integral Characteristics of Interior 
Transformations and Pseudo-harmonic Functions. III. Causal Isomorphisms in the 
Theory of Pseudo-harmonic Functions.’’ Papers I and II are to appear in the Annals 
of Mathematics. 

Morse, M., ‘“‘The Topology of Pseudo-harmonic Functions.’’ Completed but not 
yet submitted for publication. This paper obtains the basic topological relations A 
upon which Paper II is based. Paper III requires the use of infinite groups, and obtains 
relations A as a special case. Relations A can be obtained of themselves more simply 
and are so obtained in the present paper. 
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TOPOLOGICAL METHODS IN THE THEORY OF FUNCTIONS OF 

A COMPLEX VARIABLE. II. BOUNDARY VALUES AND 

INTEGRAL CHARACTERISTICS OF INTERIOR TRANSFORMA- 
TIONS AND PSEUDO-HARMONIC FUNCTIONS* 


By Marston MorRSE AND MAURICE HEINS 
INSTITUTE FOR ADVANCED STUDY 


Communicated July 25, 1945 


1. Interior Transformations.—The objective of the authors is to find 
out what properties of meromorphic functions are essentially topological 
in character, and by the introduction of appropriate topological concepts 
and methods to extend and simplify the present theory. To that.end it is 
natural to replace the class of meromorphic functions by their topological 
generalizations, the interior transformations, studied by Stoilow,! Why- 
burn? and others, and the class of harmonic functions by a similar topo- 
logical generalization termed by us pseudo-harmonic functions. The point 
of the present series of papers is by no means the generality obtained in 
this way but rather the insight and control of the theory, and the new 
problems and results which are brought out. 

Let w = F(t) ¥ const. be a transformation of a neighborhood N in the 
t-plane of a point ¢) into a neighborhood Np of wo (possibly the point at 
infinity) where F(t) is analytic at fp or has a pole at &. Let ¢ = ¢(z) map 
Z onto t and N homeomorphically onto a neighborhood N;, of 2%. Then 


w= F(¢(z)) = f(z) 


will be termed an interior transformation of N, into No. The point 2 will 
be termed a zero or pole of f(z) if to is a zero or pole, respectively, of F(¢), 
and the order of 2 as a zero or pole of f(z) will be taken as the order of to as 
a zero or pole of F(t). If the inverse of F(t) has a branch point of the mth 
order belonging to the pair (to, wo) the inverse of f(z) will be said to have a 
branch point of the mth order belonging to the pair (20, wo). 

We shall be concerned with a limited connected region G in the z-plane 
bounded by v Jordan curves 


(B) iad (Bi, vey B,). 


Without loss of generality these curves can be taken as circles. We shall 
admit transformations w = f(z) of G (the closure of G) into the extended 
w-plane such that w = f(z) is an interior transformation of some neighbor- 
hood of each point of G and is continuous on G at each point of (8). Such 
a function f(z) will be called an interior transformation of G. 

The “‘integral characteristics’’ of f(z) are as follows: (0) = the number 
of zeros of f(z) on G; n(«) = the number of poles of f(z) on G; » = the 


. 
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number of ramification points with antecedents on G of the inverse of f(z). 
These points are counted with their multiplicities. 

The order q of a sensed closed curve g in the w-plane with respect to 
w = 0 is well defined provided g does not pass through w = 0. If g is 
locally simple in the sense of Morse and Heins, I (these PROCEEDINGS, 
preceding article), the angular order p of g is well defined. A basic theorem 
follows. 

I. If w = f(z) is meromorphic on G and continuous at each point of (8), 
if the images g; of the positively sensed boundary curves 8; are locally simple 
and do not tntersect the origin, then 


n(0) + m(o) —-wh22—v+2%q — =p (1.1) 


summing over the boundary images g;. 

Examples exist in which the equality in (1.1) is excluded. Stoilow” ® has 
considered the case of one boundary curve whose image is a Jordan curve 
and has mistakenly affirmed that 


n(0) + n(o) —yw=1. 


For the Stoilow case p = gq = 1, so that the right member of (1.1) is 1 but, 
as stated, the inequality need not hold. 

II. If, however, one adds the hypothesis that the transformation w = f(z) 
1s 1-1 in some neighborhood relative to G of each boundary point of G, then 
(1.1) becomes an equality. 

The theorem with the equality has been proved by us for the case of 
interior transformations, the relation (1.1) in I only for meromorphic func- 
tion. The proof of I depends upon the lemma that the antecedents of 
points of ramification cluster at no boundary curve 6; whose image g; is 
locally simple, and this has been proved only when w = f(z) is meromorphic. 

A generalization of a theorem of Rad6* is given among various applica- 
tions. 

2. Pseudo-harmonic Functions.—If f(z) is an interior transformation of 
a neighborhood N of zo without pole at 2, then Rf(z) = U(x, y) will be said 
to be pseudo-harmonic on N. In case f(z) has a zero or pole at 2 but is not 
identically zero, then 


log| f(z)| = U(x, ») (2.1) 


will be said to be pseudo-harmonic on N except for a logarithmic pole. 
Functions U(x, y) will be admitted on G which are pseudo-harmonic on 
some neighborhood of each point of G except for logarithmic poles and which 
are continuous on G at points of (8). 

Saddle points of U(x, y) on G are similar to saddle points of harmonic 
functions. To treat the boundary values we begin with Boundary Condi- 
tion A by virtue of which the function U defined by U on (8) has at most 
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a finite number of points of extremum. It is seen that the locus U = c 
at a given level c consists of a finite number of closed arcs, or arcs which tend 
to definite points of (8). Neighboring an interior point or a boundary 
point which is not isolated at the level c, the regions “‘below c’’ alternate 
with the regions “‘above c’’ in a well-defined manner. The multiplicity yp 
of a saddle point P is one less than the number of these ‘‘sectors’’ below c 
neighboring P. The interior or saddle points of positive multiplicity to- 
gether with the points of relative minimum of U are called the critical 
points of U. 

Set M = the number of logarithmic poles of U, S = the number of 
saddle points of U on G, s = the number of boundary saddle points, 
m = the number of points of relative minimum of U, counting the saddle 
points with their multiplicities. Then 


M-S=2-v+s-—~m. (2.2) 


It is noted that U(x, y) may be without partial derivatives both on G and 
G, and that if these partial derivatives exist they may be null on (6). 
Under such conditions the classical alternatives to our evaluations, such as 
the Kronecker integral index, are undefined and therefore inapplicable. 

Boundary Conditions B.—There is a special evaluation of s — m on the 
right of (2.2) when U is of class C’ on a neighborhood of (8) and has a non- 
null gradient \ on (8). These conditions are termed Boundary Conditions 
B, and do not imply that U is differentiable in general on G. 

To proceed one covers the points of (8) at which \ is normal to (8) and 
entrant into G, by a finite set 2 of non-intersecting open arcs which exclude 
from their closures each point of (8) at which A is normal to (8) and exeunt 
from G. 

Let A(P) be the projection of \ onto the tangent to (8) at the point P 
on (8). 

The Vector Index J of (8).—Let h be any are of 2. If h is a curve 
8, an index 0 is assigned toh. If his an are of Q not a curve 6;, an index 
1 is assigned toh. To each end-point P of h at which \(P) is entrant rela- 
tive to 4 an index 1 will be assigned. To other end-points 0 will be as- 
signed. 

The vector index J; of U on B; is defined as Ey — E,, where Ep is the sum 
of the indices of end-points of arcs of 2 and £, is the sum of indices of the 
ares of 2. To (8) will be assigned the vector index 


J = J; 


A first theorem is that J is independent of the choice of the covering 2 among 
admissible coverings. Then one readily proves that 


M—-S=2-v7v4+J 
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for the most general pseudo-harmonic function satisfying Boundary Condi- 
tions B. In the special case in which Boundary Conditions A and B both 
hold, J equals the number of points of maximum of U at which is entrant 
minus the number of minimum of U® at which d is entrant. 

3. Comparison Theorems for Harmonic Functions.—We begin with the 
case where U(x, y) is harmonic on the unit circle x? + y* < 1 and continuous 
on the boundary. If u(r, 6) represents U in polar coordinates, it is found 
that when the boundary values u(1, 6) = (6) are of class C’ and p’(6@) 
satisfies a Hélder condition, then u,(1, 0) exists, and 


1 a 0 = a ° 
u,(1, ¥ ua Se #6 shoe P( | da [if b’ (80) a2 0}. 


A boundary point @ is called entrant if u,(1, 8) < 0 otherwise exeunt. The 
importance of these distinctions has already been indicated. 

A function p(@) of the above sort will be termed R-admissible if it has at 
most a finite number of relative maxima and minima occurring at isolated 
points between which p’(@) is never 0. 

If p(@) is R-admissible, has NV absolute minima and N absolute maxima 
and no other extrema, then the harmonic function with boundary values 
p(@) has just N — 1 saddle points on the région r < 1. 

In any case the number S of saddle points on r < 1 is at most N — 1. 
If p(@) is R-admissible and has 2N points of extremum and if the corre- 
sponding function u(r, 6) has S< N — 1 saddle points on r < 1, and no 
differential critical points on r = 1, it is possible to successively modify 
p(@) on arbitrarily small neighborhoods of the extremum points so that 
p(@) remains R-admissible and none of its extremum values and points are 
changed, while the number of saddle points of u (r, @) on the region r < 1 
assumes each integer from S to NV — 1, inclusive. 

If p(@) is R-admissible, it is possible to admissibly modify p(6@) on an 
arbitrarily small neighborhood of any one of its points of relative minimum 
(maximum) so that the resulting harmonic function has no critical points 
on the region r < 1. 

Various other comparison theorems are given for harmonic functions 
with and without logarithmic poles. 





* Morse, M., and Heins, M., ‘‘Topological Methods in the Theory of Functions of a 
Single Complex Variable: I. Deformation Types of Locally Simple Plane Curves. 
II. The present paper. III. Causal Isomorphisms in the Theory of Pseudo-harmonic 
Functions.’’ Papers I and II are to appear in the Annals of Mathematics. For an 
abstract of Paper I see these PROCEEDINGS, 31, 299-301 (1945). 

Morse, M., ‘‘The Topology of Pseudo-harmonic Functions.’’ Not yet submitted for 
publication. 

1 Stoilow, S., (a) ‘Du caracté-e tovologique d’un thésréme sur les fonctions méro 
morphes,” C. R. de l’Acad. des Sci. de Paris, 190, 251-253 (1930). (b) Legons sur les 
principes topologiques de la théorie des fonctions analytiques, Paris, 1938. 
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2 Whyburn, G. T., Analytic Topology, New York, 1942. 
3 Rado, T., “Zur Theorie der mehrdeutigen konformen Abbildungen.” Acta Szeged, 
1, 55-64 (1922). 


ABSOLUTE SCALAR INVARIANTS AND THE ISOMETRIC 
CORRESPONDENCE OF RIEMANN SPACES 


By T. Y. THOMAS 
DEPARTMENT OF MATHEMATICS, INDIANA UNIVERSITY 
Communicated August 13, 1945 


Let R, be a Riemann space (defined by a positive definite quadratic ' 


differential form). Denote by h,...,J, with 1 S p S ma set of absolute 
scalar differential invariants of R,. These invariants determine absolute 
scalar functions J;(x),..., I(x) of the coordinates x!,...,x” of R, with 


functional matrix 


an || anil 

‘Ox! A 

J creee teres if 

! nie miss Sipkaoe i 

Ol» oly 

(Ox! ox" 
If this matrix has rank p at all points of R,, the scalars J;, . . ., J, are func- 
tionally independent, and if the functional matrix of every q > p scalar 
invariants has rank not exceeding p the above scalars J, ... , J, (whose 


functional matrix does not vanish over R,,) are said to constitute a funda- 
mental set of scalar invariants of R,. In the following we deal with condi- 
tions for the isometric correspondence of two Riemann spaces on the 
basis of their absolute scalar invariants. The discussion is of a local nature 
so that the term Riemann space as here used is identical with the concept 
of the neighborhood of a point of a Riemann space when the im Grossen 
viewpoint is adopted. The method depends primarily on the functional 
independence of scalar invariants rather than on the order of the differen- 
tial invariants from which these scalar functions are derived. Differen- 
tiability requirements on the components g,,(x) of the fundamental tensor 
of R, are immediately obvious from the differential invariants employed in 
any instance and will not be stated explicitly. 

A space R, having a fundamental set of invariants ,..., J, with p 2 
1 will be said to be of category p. Obviously p cannot exceed n. In par- 
ticular if all scalar differential invariants of R, reduce to constants as func- 
tions of the coordinates we said that R, is of category 0. 
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If R,, is of category p, where O S p Sn, and is in isometric correspondence 
with an R,, then R,, is of category p. For if p = 0, R,, is a space of constant 
curvature. Hence R, must be of constant curvature and from this fact it 
follows readily that all scalar functions [(x) of R, are constants. Hence 
R,, is of category 0. In general the p(= 1) independent scalars J;, .. . , J, 
of R,, will be carried into p independent scalars li, .. . , I, of R, by the cor- 
respondence; hence R, cannot be of category gq < p. Similarly if there 
were q > p independent scalars in R, these wuld yield g > p independent 
scalars in R, contrary to the hypothesis that R, is of category p. 

Let li, ..., J, be a fundamental set of scalar invariants of a space R, of 
category m. For k fixed, the quantities J,, (partial differentiation with 
respect to x*) will be the components of a vector invariant of the space and 
hence the quantities 1, = g**J,,/ zg Will be absolute scalar invariants. 
Denoting by /, and I, the corresponding scalars for a space R, a necessary 
condition for the isometric correspondence of R, and R,, is the consistency 
of the scalar equations 


I(x) = (x), I(¥) = I(x) (1) 


where j,k = 1, ... ,m. Let x* = x*(x', ... , X”) be a solution of (1). 
Then the x*(x) will be differentiable functions; this follows from the 
theorem on the solution of implicit equations and the fact that the deter- 
minant |Zrce| ~ Dover R,. Now from the first set of equations (1) we have 


Se ox? 
Ljq(X) = Ijg(x) oll 
Ox’ 


From these relations and the second set of equations (1) we find 


Ox" Ox” 
pee woe noi “|=0. 
ay E ox axe F 





Since [J ta| ~ 0 it follows that the bracket expression in the above equations 
must vanish. The solution x* = «*(x) of (1) thus yields a (non-singular) 
isometric correspondence between R, and R,. Hence a space R,, of category 
n having a fundamental set of scalar invariants I,, ..., I, can be put into iso- 
metric correspondence with a space R, if, and only if, the equations (1) re- 
lating the coordinates x* and x* of these spaces, admit a solution x* = x*(x). 
If this solution exists it automatically gives the isometric correspondence. 

In what follows we limit our attention to spaces Ry. If R2 is of category 
2 the above result shows that five absolute scalars of R, suffice for the 
solution of the problem of the isometric correspondence between R; and any 
other space R,. Now suppose that R; is of category 1 and let J be the 
single independent scalar invariant of Ro, i.e., the matrix ||OJ/Ox*|| has 
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rank | over R, (in particular J may be the Gaussian curvature K of R.). 
Then form the scalars J; = a SP and Jy = i" Sz where the J,, are the 
components of the second covariant derivative of J. 

LemMaA. A space R2 of category 1 having the fundamental scalar invariant 
J admits coordinates x, y (covering any point P of Re) for which J = x and 
the line-element has the form 


ds? = p(x)dx? + q(x)dy*. (2) 


We can suppose 0//0x' ¥ 0 at any point P of Re where x', x? are the 
coordinates of a suitably chosen system. Then 7! = J(x', x?), n? = x? 
defines a non-singular transformation of a neighborhood of P. Let gg > 
6.8(n) by this transformation. Now consider the equation 

ge OF gn © 4 gir _ 9, (3) 
On® On! On? 

Since 6'! ¥ 0 this equation admits a solution v(7', yn?) defined in the neigh- 
borhood of P, such that 0v/On? ¥ 0 (see the discussion in C. Caratheodory, 
Variationsrechnung und partielle Differentialgleichungen erster Ordnung, 
1935, p. 24). Hence the transformation u = 7!, v = v(n!, n?) is non- 
singular and if 6,3, — Rag(u, v) relative to the u, v coordinates, then k!? = 0 
by (3) where the indices 1 and 2 refer to the variables u and v, respectively. 
Hence ky, = 0. Also J = wand hence J; = k!! = f(u) since J; is function- 
ally dependent on J by hypothesis. Since k'! = 1/ky it follows that 
ku = p(u). We now have 


J, = —JVigk® = —k TL, = —k'Th, — kT, (4) 


where we have used the [’’s to denote the Christoffel symbols relative to the 
u, v system. 

Ncw ae = 0. Taking a = 6 = y = | the resulting equation shows 
that I'j, is a function of u alone. Since J2 is a function of u alone, it now 
follows that the quantity k2°I'» in the right member of (4) is a function of u 
alone, i.e., 


me 1 ‘ 3 
ket, = — = o(u), (5) 


or, 


0 log ke» 
me = kud(u) = p(u)o(u) = ¥(u), 


when we substitute the expression defining the Christoffel symbol I, and 
replace the k** by their values in terms of the k,, in (5). Integrating the 











VoL. 31, 1945 MATHEMATICS: T. Y. THOMAS 309 


latter equation we see that /2 has the form ky» = g(u)r(v). Then making 
the transformation 


x=uy= fri) do 


we obtain the form of the line-element given by (2), and J = x. This 
completes the proof of the lemma. 

Now consider another space R, of category 1. A necessary condition 
for the isometric correspondence of R; and R; is that the three equations 
J=J,J: = Jvand Jz = Jz admit a solution ¥* = X%*(x!, x?) where J, J; 
and J» denote the scalars in Re which correspond to the scalars J, J; and 
J2in Re. We assume the existence of this solution. Since J; ¥ 0 in R, 
it follows that J: ~ 0 in R,. Hence J is a fundamental scalar invariant 
for R,. Now suppose P > P by the correspondence x* = X%(x', x”). In- 
troduce the coordinates x, y in the neighborhood of the point P as given by 
the above lemma for the space R, and the corresponding coérdinates x, y 
in the neighborhood of P for the space R, relative to which we have 


J = &, ds* = p(z)dx* + G(z)dy?. 


The above solution x* = x*(x', x?) implies the existence of a solution ¥ = 
x(x, y), y = (x, y) of the equations J = J, J; = Ji and J; = Jz relative to 
these new codrdinates. But from J = J we have ¥ = x. Also the rela- 
tion Ji = J; yields g"\(x) = g!(x) from which we deduce gu(X) = gu(x) 
or p(x) = p(x). Then from Jz = Je we find 


eTn + gle = g'Th + eT ns. (6) 


Since g" depends on x alone it follows that Ti: depends on x alone. Simi- 
larly for Tj, Also from the relations g'! = g", g, = gy and ¥ = x we 
obtain I'j,(x) = Ti,(x). Hence (6) reduces to 2"T% = g”T. On simpli- 
fying, this latter relation becomes 


Then using the fact that ¥ = x and integrating we have g = C’q where C is 
a constant. 

From the relations p(%) = p(x), g(%) = C’q(x) and ¥ = x we see im- 
mediately that the transformation x = x, y = Cy + d where d is a suitable 
constant gives an isometric correspondence between the neighborhoods of 
Pand P. Understanding that the term Riemann space is identical with a 
suitably restricted neighborhood in accordance with the local viewpoint 
the result obtained can now be stated as follows. A Riemann space R: of 
category 1 with fundamental scalar invariant J can be put in isometric corre- 
spondence with a Riemann space R; if, and only if, Re is of category 1 and 
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the three scalar equations J =J,J: = Jiand Jz = Jz admit a solution x* = 
x*(x!, x*) relating the coérdinates x* of Rz to the codrdinates x* of Ro. 

Since a Riemann space R; of category 0 has its Gaussian curvature K = 
const. we can immediately state the following result. A Riemann space 
R, of category 0 can be put in isometric correspondence with a Riemann space 
R, if, and only if, Re is of category Oand K = K where K and K are the Gaus- 
sian curvatures of R,and Ro, respectively. 1t has already been observed that 
the categories of two Riemann spaces in isometric correspondence must be 
the same. The remainder of the proof consists in showing that two Rie- 
mann spaces of the same constant curvature can be put in isometric corre- 
spondence and this is a well-known result. 

We leave open the problem of the extension of these results to Riemann 
spaces R, where » > 2. Anvther problem, worthy of serious study, con- 
cerns the investigation of the isometric correspondence of Riemann spaces 
in the large (for example, open and compact spaces) on the basis of their 
absolute scalar invariants. 


BERNOULLI’S NUMBERS AND CERTAIN ARITHMETIC 
QUOTIENT FUNCTIONS 


By H. S. VANDIVER 
DEPARTMENT OF PURE MATHEMATICS, UNIVERSITY OF TEXAS 
Communicated July 6, 1945 


In another paper! the writer established a congruence involving Ber- 
noulli numbers, which can be extended, as we shall now show, to the 
generalized Bernoulli numbers of the first order. Such a number, b,(g, 4), 
is defined as the number obtained by expanding (gb + h)” in full by the 
binomial theorem and substituting , for b* where b, is obtained by the use 
of the formula (6 + 1)* = },,7 > 1, and the left-hand member is expanded 
as before and b, is substituted for b*. Let pb be an odd prime. Then 
pb — 1 has the form mc, with c even. Now consider the expression 


X(1 — (dp + 1)?-?)? = N (mod fp), (1) 


where p ranges over all the values such that p° = 1 (mod p?). Any term 
in this expression reduces to unity, modulo p’, if p is such that 


dp + 1=0 (mod p). (2) 


For each p that does not satisfy (2), the corresponding term of (1) is di- 
visible by p*. In the relation x° = 1 (mod p”), the c incongruent solutions 
modulo p? are also incongruent modulo p. Hence, N in (1) is the number 
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of incongruent solutions p, modulo p, in (2) and hence equals 1 or 0, ac- 
cording as d° = 1 (mod p) or d° # 1 (mod p). Expansion of (1) gives 


Nac — 2% do d® ‘a ty +o ae (70 | z ") 
k=0 Re 


modulo p?, and this may be written 
on ee . ke eee . ke 2(p — 1) 
N ae Sat todd ( Mie 


k=1 
. m 2(p 7 7) 
(m+k)e 
cx 4 be + kc} 
Set d = ga + hand multiply through by (ga + h)", with m # 0 (mod c). 


Using the value of N noted above, we have, after letting a range over the 
integers 0, 1, 2, ...,  — 1 and adding, 


u (gr + h)* = —2c a Sre+n(Z, » i a ') é 
: E. Steals n(2, 2 = te Pa Smetretn(8s h) oh 1 ee (3) 
modulo p?, where 


p-l1 
S(g, h) = » (ga + h)' 


and r ranges over the integers in the set 0, 1, 2..., p — 1 with (gr + h)° 
=1(mod p). It can be shown, as follows, that for 7 even and p > 3 


Si(g, h) = pbi(g, h) (mod p?). (4) 
As noted in. another? (relation (6), page 576 of the first reference) paper, 





_ifitl 
i + DSde,¥) =O (E41) dee wpertgte 
which gives 
4 am 
Sieh) = D s(¢ Voena(es Bp PO pig. (5) 
i 
For gq > 1 


p> +2) 214+2(¢-)=%¢™W-lLegqtl. 


Hence each term of (5), beyond the second, is divisible by p?, since the factor 
p cannot occur in the denominator of any fractional part oftener than in 
the numerator. But the second term will also be divisible by p, for 7 is 








312 MATHEMATICS: H. S. VANDIVER Proc. N. A. S. 


even and b,_1(g, ) is an integer.? This leaves (4). For 7 to be even in (4), 
n must be even, and the limitation p > 3 will be met; because if p = 3, 
then, as p — 1 = cm,c = 2, butm #0 (modc). Using (4) and dividing 
(3) through by p, we obtain, modulo p, with c and even, 


> (gr + h)” pet 
— (2(b — 1) ™ (Ap — 1) 
Cc 2» ( hes ) Dec +nlZ, h) + ch ‘oe + sa Dme+ke+n(2> h). (6) 


Reducing the binomial coefficients, modulo p, as on page 58 of a previous 
paper,! we obtain modulo p, 


as eee = —2¢ (—1) sad EN 2 h) vA 


c > (—1)¥(Re + 1)Die nlg, 2) +6 (— 1) heby rt ect nl 2). (7) 
k=1 k=] 
Now use the relation® 


bi+- i(g, h) _ ws h) 
s+tp-1 


where g # 0 (mod #) and for s # 0 (mod wo — 1), which latter limitation is 
here complied with since » # 0 (mod c). We may then write, modulo p, 
with g # 0 (mod 9), 


Ler +h” im 
ida cain =C pm (ke + 1) dic4n(2, h) ¥ 


k=1 


(mod p) 





: 


m ketn—1 = bic+-n(g, h) h) 
ans Sw Terreeganed c+n\s>» —2 c+n\S . 
ch ke "aa bicyn(Z, 1) oD ductal h) = —cn ote ag ge 


This gives 
THEoREM I. If p is an odd prime with p — 1 = mc, c even, and g # O 
(mod p), then 


> (gr + h)” m—1 
r se Dxc+-n(Q, | h) 
ge cn » ee a (mod p) (8) 





for n even and not divisible by c, where the summation on the left extends over 
all integers r in the set 0,1, 2, ..., p — 1, such that (gr + h)’ =1 oe p). 
For the special case g = I, h = 0 this.is Theorem I of the previous’ paper, 
and we shall now proceed to find another congruence involving M(c, n, p), 
that symbol being defined as the left-hand member of (8) for this special 
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case. In another‘ paper the writer defined a repetitive set modulo / as a 
set of integers 
he ee (9) 
where / is any integer and r; < /, such that there exists an r # 1 with 
Fg Wag 0s ig 1% 


being congruent in some order to (9) modulo /. Hence, for / = p, the 
integers r in the set 1, 2, ..., ’ — 1, which satisfy r° = 1 (mod fp), forma 
repetitive set modulo p with any such r # 1 as multiplier. As the writer 
proved, the set of integers 


Yh Ts, gm 1, 8 


ey 


where 


yZ=- 5 (mod 1); 0OS%<1r 
is a permutation of (9). Hence 


) (1e ts ms re) = Dri + np D yer 3-* (mod p’) 


a 


or 
(r® — 1)M(c, n, p) =n D> yer *~* (mod fp). (10) 
If 
a 

Y= > (mod r),0 S¥< 7, 
then for the special case c = p — 1 (10) gives 

$=—1 = 

n > yea"! = (r" — 1) > (mod p), 

a-—l 

and it follows that 


Sygate Gein), (11) 


1 
a=1 n 
for n # 0 (mod p — 1). 


Let n = cs + e,0< e<c, where s ranges over the values 0, 1,...,m — 1, 
and let r belong to the exponent ¢c modulo p. (11) gives on addition 


m—1 
= ( 5 Le Te uae Dr ma 
s=0 CS + é sa s 4a 
(12) 
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where a; ranges over the a’s satisfying a© = 1 (mod p) and a over thea’s 
satisfying a # 1 (mod p). 


We have 
Dd ate? = af? (1 +aj tay + ....+a{"—) =a{-' m (mod ), 


$s 
and 


Dag =a (1 +a +axy+.... +a") = 
eee 
a,” os = 0 (mod p) 


ay — 


since a3 # 1. Hence, the right hand member of (12) reduces modulo p to 


1 


myo 0i 
a, 
Multiplication by —ce gives 


m—1 b 

—ce(r’ — v( > Ste | D eva ai*. (13) 
s=0 cS +e a 

But we also have (10). Since r belongs to the exponent ¢ modulo p, and 

0 <e<c, and 


Ost oma a bs (mod p) 
s+p-! s 


for s 0 (mod p — 1), comparison of (13) and (10) also yields the above 
noted special case of Theorem I for g = 1, h = 0. 

It is not clear how the argument used in this second proof of Theorem I 
of the former! paper can be extended to give a proof of Theorem I of the 
present paper, since if we attempt this extension it seems necessary to 
employ a relation in another paper‘ (last separate congruence on page 123) 
and that relation was proved under a number of restrictions which would 
not apply in the statement of our general theorem. 


1 These PROCEEDINGS, 31, 55-60 (1945), Theorem I. On page 55 of this paper on the 
second line above (4), insert the factor (—1)" in the right-hand member. The relation 
(4) is subject to the condition  # ku. On page 56, second line, read ‘‘are at most only 
(pb — 3)/2” in lieu of “‘are only p — 1.” In the first line above relation (6), insert 0 after 
the second congruence sign. In the eighth line above relation (7), read a < in lieu of 


a <(p— 1). On page 59, fourth line above relation (14), read Zr” in lieu of Zr. In 
r r 


connection with the last formula on page 60, add the condition (n(m + c), p) = 1. 

2 These numbers were considered by the writer in Duke Math. Jour., 8, 575-584 
(1941), and also in Trans. Am. Math. Soc., 51, 512-519 (1942). 

3 These PROCEEDINGS, 28, 326 (1942), relation (8) for 7 = 1. 

* Bull. Am. Math. Soc., 46, 122 (1940). 











VoL. 31, 1945 MATHEMATICS: A. D. MICHAL 315 


THE TOTAL, DIFFERENTIAL EQUATION FOR THE 
EXPONENTIAL FUNCTION IN. NON-COMMUTATIVE NORMED 
‘LINEAR RINGS 


By ARISTOTLE D. MICHAL 
CALIFORNIA INSTITUTE OF TECHNOLOGY 


Communicated June 30, 1945 


A normed linear ring! is a normed linear space (say with real multipliers) 
over which there exists an associative (but not necessarily commutative) 
function AB (additive and continuous in each variable) with values in 
the space. We shall assume the existence of a unit element I and assume 
that the space is a complete normed linear space (Banach space). We need 
not normalize the modulus of AB and take it as the number 1, nor need 
we normalize the norm of the unit element J and take it as the number 1. 
Examples of normed linear rings are sets of square matrices, quaternions 
and more generally any (real) associative algebra topologized with a norm 
topology in a variety of possible ways. Important infinite dimensional 
examples abound in the theory of integral‘equations and functionals. 

The exponential function e4 in a normed linear ring is defined by 


A? 
I+ At a + oes + ass 


Since ||AB|] < m||A|| ||B||, where m is the modulus of the bilinear function 
AB, it follows that the exponential function in a normed linear ring is an 
entire analytic function.2 The object of this paper is to characterize the 
exponential function in a normed linear ring by a total differential system 
in Fréchet differentials. More specifically, this characterization is em- 
bodied in the following theorem. 

TuHEorREM. The non-linear total differential system in Fréchet differentials 
of functions 2(A) with arguments and values in a normed linear ring N 


i(A) = Ji'e((h— DA)dda(eA)ae A 
2(0) = J, the unit of N, 


has a unique entire analytic solution given by the exponential function 2(A) = 


e4 in N. 
Let 


2(A) = >52,(A)(z,(A), a homogeneous polynomial of degree 7) 
i=0 


be an entire analytic solution of (1). Clearly 29(4) = J and by calculation 
2,(A) = A, (A) = A?/2!. In this and in what follows, a theorem by the 
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author* on term by term Fréchet differentiation of a power series in normed 
linear spaces is used. The following recurrence formula holds for n > 2 
and for all A, 64 e N: 


(n+ 1)Py4:(A,A, ...,4,4) = Jf, '{042,(EA) + 2,((1 — £)A)8A + 
Yei((1 — 8)A)6Az(€A)}dé, (2) 


t+j=n 
4,721 


where P,,4,(A1, As, ..., An+1) is the polar of the polynomial z,,4,(A4). The 
truth of the theorem follows after some calculation with the aid of the 
recurrence relation (2). 

The ideas that led the author to the differential system (1) for the par- 
ticular case in which the normed linear ring N is a matric ring are discussed 
in the author’s paper, ‘Differential Equations in Fréchet Differentials 
Occurring in Integral Equations,” these PROCEEDINGS, 31, 252-258 
(1945). 

The characterization of other elementary functions in a normed linear 
ring NV, such as sin A and cos A, can also be effected by total differential 
systems in Fréchet differentials. Similar results can also be given for 
normed linear rings without a unit. One can, for example, study the 
properties of the function defined by the expansion 


 . A 
Ato tat + 


It should be pointed out that the differential system (1) is not of the 
type for which Michal and Elconin have given existence and uniqueness 
theorems.‘ 


1 A search through the literature seems to reveal that the first systematic study of 
normed linear rings (non-Hilbert variety) was made by Martin, R. S., and Michal, A. 
D., in ‘Some Expansions in Vector Space,” Jour. Math. Pures et A ppl., 13, 69-91 (1934). 
This paper was presented before the Los Angeles 1932 national meeting of the American 
Mathematical Society. The author has made an extensive use of normed linear rings 
in his studies on general differential geometry. See Michal, A. D., “General Differential 
Geometries and Related Topics,’’ Bull. Amer. Math. Soc., 45, 529-563 (1939), where 
earlier references are included. See also Michal, A. D., and Mewborn, A. B., “Abstract 
Flat Projective Differential Geometry,’’ Acta Mathematica, 72, 259-281 (1940). For 
problems in mathematical analysis, besides those treated in the Michal-Martin paper, 
the reader is referred to Michal, A. D., and Hyers, D. H., ‘“‘Second Order Differential 
Equations with Two Point Boundary Conditions in General Analysis,”” Amer. Jour. 
Math., 58, 646-660 (1936); and to Michal, A. D., and Elconin, V., ‘‘Completely In- 
tegrable Differential Equations in Abstract Spaces,’ Acta Mathematica, 68, 71-107 
(1937). The last paper also considers briefly some questions on ideals in normed linear 
rings. In recent years many important coritributions have been made to the theory 
of normed linear rings by many authors in Russia as well as in America. The Russian 
mathematician I. Gelfand is especially to be mentioned in this connection. See his 
paper in Rec. Math. (Mat. Sbornik), 51, 3-24 (1941). For commutative normed linear 
rings, see Lorch, E. R., Bull. Amer. Math. Soc., 50, 447-463 (1944). 
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2 Cf. Michal-Martin paper, loc. cit. 

3 Michal, A. D., A Theorem on the Fréchet Differentials of Power Series in Normed Linear 
Spaces (unpublished). 

4 Cf. Michal-Elconin paper, loc. cit. ‘ 


UNIVERSAL RATIONAL FUNCTIONS 
By E. T. BELL 


DEPARTMENT OF MATHEMATICS, CALIFORNIA INSTITUTE OF TECHNOLOGY 
Communicated July 17, 1945 


1. Positive in all that follows shall mean greater than zero. The coeffi- 
cients in all rational functions considered are integers (actually all positive 
in the examples given). 

According to a customary definition, a form F is universal if, for integer 
values (positive, zero or negative) of all the variables in F, F represents all 
integers, positive, zero or negative. Another type of universality, as in 
Waring’s problem, requires that a form represent all positive or zero 
integers for positive or zero values of the variables. Thus the sum of four 
squares is universal in this sense, but is not universal if zero values of the 
variables are excluded. A further type restricts the universality to some 
proper subset C of all the integers, for example, cubes. 

We shall call the rational function F/G, where F, G are forms with 
positive integer coefficients, and F/G is in its lowest terms, universal, if, 
for positive integer values of all the variables in F, G, F/G represents all 
positive integers. If F/G for positive integer values of the variables 
represents all integers in the subset C, F/G will be called universal for C. 
The universal rational functions considered here also represent all negative 
integers when negative integer values of the variables are admitted. These 
particular universal rational functions therefore represent all integers 
except zero for integer values of the variables different from zero. A point 
of interest, however, is that positive integer values of the variables suffice 
for the representation of all positive integers. 

2. An example of a universal rational function in four variables x, y, 
Z, wis 

x2(x"—MyMt1 4 gh—lyn+1) 

gg ts oe w"—ignt+l 

where » is an arbitrary constant positive integer; 
eetlghtigtt: + A, his, stale 
emLyn— ly sl gm —layt—Iyt—1 


in six variables x, y, 2, u, v, W, is universal for squares, where m, n, s are 
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arbitrary constant positive integers. In 2s variables x, ..., %;, Vi, ..-» 
ys, with the m,, n, arbitrary constant positive integers and 


s s 
F(x, s++y Mey Vy -- Ys) =I eget! + Il ett, 
i=l ; 


i=1 


s s 
Gn . »<9 Bee Pas + +9 J) SU a + Hy s—', 
i=1 


i=1 
F/G is universal for squares. With s > 1 and 
s—1 s—1 


F =x, 0 x*! + y, I y%*, 


t=1 i=1 


s—1 s—1 


G =y, I x%-! + x, I y,"-, 


i=l i=l 
F/G is universal. For even integers, 
xr? + ys? + 274? + Buu? 
xr + ys + 2t 


is universal. For fourth powers, (x° + u*v*w)/(x + w) is universal. 
Illustrative of an infinity of universal rational functions F/G in which 
each of F, G is a binomial, x?w?(y? + zw)/(x> + w*) is universal. 

3. The foregoing examples suggest that universal rational functions 
are much more easily found or constructed than are universal forms. That 
this is so, is evident from the remark that an elementary monomial repre- 
sents all positive integers for positive integer values of the variables: the 





monomial x;'. . .x,%, in the s independent variables x, ..., x;, is elementary 
if at least one of ihe positive integers a;, ...,@,is 1. Without loss of gener- 
ality, we shall assume a; < d2 < ... <@,; (Gi, @, ..., ds) is then the index 


of the monomial. 

If F/G is universal, F = nG, where n is an arbitrary positive integer, has 
a solution for positive integer values of all the variables y,, ..., y, in F, G. 
Suitable F, G are found by representing m as an elementary monomial, say 
n = my”. ..ns%, of index (a1, ...,@;). For a given index, the total number 
of such representations is readily determined as a combinatorial function 
of the exponents in the prime decomposition of m. The representation of 
n in F/G is then of the type 


y= M,(ni, eens Ns), 4 <= i; oe 2 


where the VM, are monomials in 1, ..., 3. 
Write F/G = R(y, ..., ¥,). From what precedes, it follows that 
R(%, ..., 3)R(wy, ..., W:) = Rew, ..., BW) 


has an infinity of solutions in monomials 2, ..., 2, W:, ..., W: For if 
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nN = Mm%...ns8, mM = m"...m,% are representations of the arbitrary 
positive integers m, m in a monomial of index (a, ..., @;), their product is 
represented as (mym)"...(m,m,)% in a monomial of the same index. If 
2, = Mi(m, ..., ms), Ww = M,(m, ..., ms) represent n, m, respectively, 
in R, 


zw, = M(nym, ..., nsm,) 


for the representation of mm in R. 


It follows also that if R(x, ..., x;) is a universal rational function, and 
Ry, ..., Ry are any universal rational functions in any variables, then 
R(Ri, ..., R;) is a universal rational function. For each variable denotes 


a positive integer, and hence each is representable in any universal rational 
function. 

Appropriate F, G may be constructed indefinitely from the disjunction 
(logical product) and conjunction (logical sum) of simple multiplicative 
diophantine equations. If the totally distributed form of the disjunction 
or conjunction can be written as F = MG (or asG = MF), where M is an 
elementary monomial, and M, G (or M, F) have no common factor, then 
F/G (or G/F) is universal. The solution of F/G = n (or of G/F = n), 
where 1 is an arbitrary positive integer, is then obtained by expressing n 
inthe form M. If F = MG (or G = MF) was obtained by disjunction, the 
values of the variables are found by solving any one of the equations in 
the disjunction; if the equation was obtained by conjunction, the equations 
are solved simultaneously. Since in either case the solutions are the 
general solutions, given in a necessary and sufficient number of independent 
integer parameters, the corresponding M represent all positive integers. 
The first example in the preceding section was constructed from the dis- 
junction of xy" = zw", yz" = wx"; the fifth, from the conjunction of xr = 
ys = zt = wu. Both may be considerably generalized. 

















